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I .  INTRODUCTION 


When  a  ship  travels  with  constant  velocity  along  the  surface  of 
a  liquid,  it  creates  behind  it  a  wake  which  is  called  a  "ship  wave" 
pattern.  A  similar  pattern  is  produced  by  a  submerged  object  moving 
parallel  to  the  surface.  The  usual  analyses  of  such  patterns  apply 
to  liquids  of  uniform  density  in  which  only  one  type  of  propagating 
wave,  called  a  surface  wave,  is  possible.  We  shall  consider  ship 
wave  patterns  in  horizontally  stratified  liquids  in  which  one  or  more 
propagating  internal  waves  exist  in  addition  to  the  surface  wave. 

Keller  and  Levy  (Ref.  1)  have  shown  that  in  any  such  liquid  the  ship 
wave  pattern  is  a  superposition  of  separate  patterns,  one  for  each 
propagating  internal  or  surface  wave.  They  have  also  obtained  formulas 
for  the  wave  height  and  particle  velocity  as  functions  of  position 
throughout  the  pattern.  From  these  formulas  one  can  see  that  for  a 
submerged  object  the  patterns  corresponding  to  some  of  the  internal 
waves  can  have  larger  amplitudes  than  that  corresponding  to  the  sur¬ 
face  wave.  Therefore  we  shall  examine  the  internal  wave  patterns  in 
detail  for  a  simplified,  tut  realistic  density  profile  in  which  in¬ 
finitely  many  propagating  internal  waves  occur.  Previously  Hudimac 
(Ref.  2)  studied  the  special  case  of  a  two-layer  fluid  in  which  just 
one  propagating  internal  wave  exists. 
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II. 


WAKE  GEOMETRY 


To  describe  the  wake  of  a  horizontally  moving  object,  we  replace 
the  object  by  a  point  which  we  call  the  source.  (  We  introduce  cartesian 
coordinates  in  the  horizontal  x,z  plane  containing  the  source,  with 
the  x-axis  along  the  path  of  the  source  and  the  origin  at  the  position 
of  the  source  at  time  t  =  0.  If  the  speed  of  the  source  is  -v  then 
the  coordinates  x^t7),  zQ(t')  of  the  source  at  time  t'  are 

xQ(t')  =  -vt ',  zQ(t ' )  =  0  (2.1) 

We  wish  to  determine  the  wake  corresponding  to  waves  of  a  partic¬ 
ular  type  emitted  by  the  source,  i„e.,  to  the  surface  wave  or  to  the 
n-th  internal  wave.  We  suppose  that  the  source  emits  waves  of  this 
type  with  all  frequencies  u>  and  that  the  wave  has  a  definite  propaga¬ 
tion  constant  or  wave  number  k.  It  is  convenient  to  express  w  as  a 
function  of  k, 


u>  =  iu(k)  (2.2) 

The  functional  relation  (2.2)  is  determined  by  the  density  profile, 
and  will  be  considered  later. 

Let  us  consider  the  phase  -0(x,z,k,T)  at  the  point  x,z  at  time 
t  =  0  of  the  wave  of  wave  number  k  emitted  by  the  source  at  time  -t, 
t  s  o.  If  the  wave  is  emitted  at  phase  zero  then 

-0(x,z,  k,T)  =  kr  -  u)t  (2.3) 
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Here  r  is  defined  by 


r  =  {[x-xo(-t)]2  +  z2}^  (2.4) 

We  seek  those  values  of  k  and  t  for  which  0  is  stationary.  This 

requirement  yields  the  the  two  conditions* 

0  =  -0^  =  r  -  u>kT  (2.5) 

0  =  -0  =  kr  -  u)  (2.6) 

From  these  equations  we  find  that 

r/r  =  c  =  u),  (2.7) 

'  g  k 

v(x  -x)/r  =  c  =  u>/k.  (2.8) 

Here  we  have  introduced  the  group  velocity  cg  and  the  phase  velocity 

c  defined  by  the  last  equalities  in  (2.7)  and  (2.8).  Equation  (2.7) 

shows  that  the  wave  from  the  source  xq(-t),0  travels  to  x,z  at  the 

group  velocity  c  .  Equation  (2.8)  shows  that  the  trace  on  the  x-axis, 
9 

of  the  straight  line  perpendicular  to  the  ray  from  xq(-t),o  to  x,z, 
travels  with  the  source  velocity  -v. 

The  two  equations  (2.7),  (2.8)  determine  the  values  of  k  and  t 

which  make  0  stationary.  When  these  values  are  used  in  (2.3),  (2.3) 

will  yield  the  stationary  va.1ue  of  the  phase  at  each  point  x,z.  These 

results  are  just  Equations  (11.5)  and  (11.6)  of  Ref.  1,  which  we  have 

redr  rived  in  a  simpler  way.  We  now  use  (2.7)  to  write  t  =  r/c  and 

y 

(2.8)  to  write  m  =  kc.  Then  we  can  rewrite  (2.3)  as 

-0  =  kr(l  -  c/c  )  (2.9) 


*Letter  subscripts  denote  partial  differentiation. 


5 


Next  (2.8)  and  (2.4)  yield,  if  c  £  v, 


2  2  "k 

xq~x  =  rc/v,  z  =  r(l  -c  /v)2 
From  (2.1)  and  (2.7), 


(2.10) 


x 


o 


w/cg 


(2.11) 


Let  us  now  eliminate  xq  from  (2.10)  by  means  of  (2.11)  and  then 
eliminate  r  by  means  of  (2.9).  Thus  we  obtain 


x  =  (0v/k) 


~1  -  (ccq/v) 


c  -  c 


1  Ml  -  (c/v)2' 

J>  z  =  L . 9  ”"c~-"c - . 


(2.12) 


Thus  if  c  £  v,  (2.12)  is  the  parametric  equation  of  the  wavefront  0  = 
constant,  where  0  is  the  stationary  value  of  the  phase  and  k  is  the 
parameter. 


It  is  clear  from  (3.4)  -and  (3.5)  that  for  long  waves,  for  which  k  is 
small,  both  x  and  z  are  large. 


To  eliminate  k  we  solve  (.3.4)  for  k  and  substitute  into  (3.5), 
obtaining 

•*  2 

U-jX  2(0W2VJx)^ 


u>. 


z  = 


,~~2  2Ti  “2  2 

(v  ~uJj)  v  -  u)j 


+  .  • 


(3.6) 


If  iu(k)  is  not  analytic  around  k  =  0,  (3.1)  is  not  valid  and 
therefore  (3.6)  does  not  apply.  This  is  the  case  for  ordinary  surface 
waves  in  water  of  constant  density  and  infinite  depth,  since  for  them 
u>(k)  =  (gk)^.  Then  (2.12)  becomes 


v  -  20v 

x  —  —  p 

(gk)* 


{''Ik)’  2^{x~k) 


(3.7) 


2 

From  (3.7)  we  see  that  k  must  be  restricted  to  the  range  k  >  g/v  in 
order  that  z  be  real,  so  a  smali-k  expansion  is  not  applicable  in  this 
case. 


8 


IV.  SHORT  WAVES  AND  THE  NEAR  WAKE 

For  short  waves  or  large  k,  we  assume  that  tu(k)  has  the  asympto¬ 
tic  expansion 


io(k)  =  N  -  N2k"2  -  N4k-4  -  ... 

(4.1) 

Then 

C  =  Nk"1  -  N2k“3  -  N4k“5  +  ... 

(4.2) 

cg  =  2N2k“3  +  4N4k“5  +  ...  . 

(4.3) 

Upon  using  (4.2)  and  (4.3)  in  (2.12)  and  (2.13)  we  obtain 

v  — 

(0v/k) 

"l- v"  2  ( Nk_1-N2k“ 2  -  . . . )  (  2N.2k"  3  + 

...)] 

A  — 

Nk"1  -  3N2k-3  -  ... 

J 

- 

(0v/N) 

L  3N2  v-2  1 

1  +  N  k  +  •  ’  j  • 

(4.4) 

z  =  (0/k) 

'(2N2k_ 

3+4N4k"5  +  ...)  (l-v"2N2k"2+2v'2NN2k“4  +  ...)?i*| 

- 

Nk-1  -  3N2k-3  -  ... 

J 

=  (20N2/Nk5) 

-  /2N,  n2  3N2\ 

vir'i7+_r/  +  --J 

(4.5) 

Solving  (4.4) 

for  k 

and  substituting  the  result  into 

(4.5)  yields 

*  =  -^^<x-r>3/2+- 

(27v3N20)^ 

(4,6) 
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From  (4.6)  we  see  that  each  wavefront  0  =  constant  has  a  cusp  at  x  = 
0v/N  on  the  path  z  =  0. 

In  the  case  of  ordinary  suz‘face  waves  in  water  of  constant  density 
and  finite  or  infinite  depth,  (4.1)  does  not  hold  so  neither  does  (4.6). 
For  infinite  depth  (3.7)  yields  for  k  large, 

z  =  ±  gx^/4v^0  +  ...  (4.7) 

Then  all  wavefronts  enter  the  origin  x  =  0  on  the  path  z  =  0.  The 
result  (4.7)  also  holds  for  the  finite  depth  case  when  the  density  is 


constant. 


V.  EXAMPLE 


A.  THE  DISPERSION  EQUATION 

Let  v(y)  be  the  y  dependent  factor  of  the  vertical  component  of 
fluid  velocity  in  a  time  harmonic  wave  of  angular  frequency  uj  and 
wavenumber  k  in  a  fluid  of  density  pQ(y)  and  depth  h.  Then  v(y) 
satisfies  the  equations 


(Ref. 

fined 


1,  Eqs.  5.14  - 
by 


+  k2((u“2N2  - 

1)  V  =  0,  O^yS-h 

(5.1) 

vy(0)  =  kV 

2gv(0) 

(5.2) 

v(-h)  =  0. 

(5.3) 

It  It  IT 

5.18).  Here  N  (y)  is  the  Vaisala  frequency  de- 


N2  =  C-g)(p0)y/p0  (5.4) 

2 

This  problem  has  nontrivial  solutions  only  if  k  is  an  eigenvalue. 

We  shall  take  for  p()(y)  the  following  function 

oQ(y)  =  ,  0  £  y  £  y2 

=  Pl  exp[(N2/g)(y^-y)]  ,  y±  *  y  ^  y2  (5.5) 

=  P2  H  pi  exp[(N2/g)(y1-y2)]  ,  y2  ^  y  ^  -h 
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o 

The  layer  between  and  y2  is  the  thermocline,  within  which  N  is 

constant,  and  outside  it  =  0.  Now  the  coefficients  in  (S.l)  are 

piece-wise  constant  so  (5.1)  can  be  solved  explicitly.  If  we  ignore 
-12 

the  term  g  N  v  in  (5.1),  the  solution  is  simply 


v(y)  =  sinh  (ky  +  a)  0  ^  y  ^  y1  (5.6) 

vfy)  «  C  cos  kvy  +  D  sin  kvy  y^  s:  y  ^  y2  (5.7) 

v(y)  =  B  sinh  k(y  +  h)  y  '<£.  y  s  _h  (5.8) 

Here  v  is  defined  by 

v2  =  (N2/u>2)  -  1  (5.9) 

Condition  (5.3)  is  satisfied  by  (5.8),  while  (5.2)  yields 

tanh  a  =  w2/kg  (5.10) 


At  y^  and  y0  both  v  and  must  be  continuous.  This  requirement 
yields  the  four  conditions 


C  cos  kvy^  +  D  sin  kvy^  =  sinh  (ky^  +  a)  (5.11) 
C  cos  kvy2  +  D  sin  kvy,,  =  B  sinh  k(y2  +  h)  (5.12) 
-v  C  sin  kvy^  +  vD  cos  kvy^  =  cosh  (ky^  +  a)  (5.13) 
-v  C  sin  kvy2  +  vD  cos  kvy2  =  B  cosh  k(y2  +  h)  (5.14) 


We  now  combine  (5.11)  and  (5.13)  to  obtain 
[C  cos  kvy^+D  sin  kvy1]cosh(ky^+Qr) 

•j-v[C  sin  kvy^-D  cos  kvy^]sinh(ky^+a)  =  0. 
Similarly  we  get  from  (5.12)  and  (5.14) 
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(5.15) 


[C.  cos  kvy2+D  sin  kvy23cosh  k(y2+n) 

+v[C  sin  kvy2-D  cos  kvy23sinh  k(y2+h)  =  0  (5.16) 

In  order  for  (5.15)  and  (5.16)  to  have  nontrivial  solutions  for  C  and 
D,  the  determinant  of  the  coefficient  matrix  must  vanish.  This  yields 
the  dispersion  equation 

tan  kv(y1-y2)[l+v2tanh  (ky-j+o?)  tanh  k(y2+h)] 

;=  v[tanh  (ky1+a)“tanh  k(y2+h)]  (5.17) 


B.  FAR  WAKE 


Let  us  examine  (5.17)  for  k  small,  tentatively  asstuning  that 

u)  ~  kc(0)  as  k  tends  to  zero.  Then  (5.9)  yields  v  ~  N/kc(0),  (5.10) 
2 

yields  a  ~  kc  /g  end  (5.17)  becomes  at  k  =  0 


tan  [(N/c)(y1-y2)3  = 


Nc“  \  yx-y2-h+ ( c  2/g ) ] 

1  +  N2c"^[y1+(c2/g)3(y2+h) 


(5.18) 


This  xs  an  equation  for  c(0)  which  has  infinitely  many  solutions  which 
we  shall  call  cn(°)j  n  ~  0,1,2,...  .  To  describe  them  we  write 

*»«>  ■  nrh;  »  (s.h> 


Here  we  have  introduced  the  themocline  thickness  s  defined  by 


s  =  Yi  -  y. 


(5.20) 


Then  (5.18)  becomes  the  following  transcendental  equation  for  a  : 

*1 


5[<s-h)<nm.an)  +  gj— 7! 


ta"  a“  =  i  r~.  ~  *  '  V2  <  an  <  n/2  (5.21) 


s^+[y^(nTt+an)  +N  s  /g3(y2+h) 
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For  n  large  (5.21)  yields 


s(h-s) 


n  nny1(h+y2J 


(5.22) 


If  h  »  s  and  h  »  lygU  (5.22)  becomes 

s  _  , 

a  ~  -  -  ,  n  »  1 

n  nny1  3 

For  n  =  0,  (5.21)  becomes 


tan  a„  = 
o 


sC(s-h)ao  +  §~-l 
2  .  N2s2 


s"  +  Cy^  +  Sf)(y2^) 


(5.25) 


(5.24) 


If  | y2 |  «  h  and  N2s/g  «  1,  we  can  replace  tan  by  a^  in  (5.24). 
The  resulting  biquadratic  equation  has  as  its  two  positive  solutions 


Ns 

°  *  (gh)% 

(5.25) 

a+  ~  (-2-)^ 

0  v-y1/ 

(5.26) 

Let  us  now  use  the  results  (5.23),  (5.25),  and 

(5.26)  in  (5.1?) 

and  introduce  the  effective  gravity  g'  defined  by 

g'  =  N2s  „  “2  '  Pl  g 

Pi 

(5.27) 

Then  (5.19)  yields 

c“(0)  -  (gh)li 

(5.28) 

Cq(o)  ~  (g'lyj )* 

(5.29) 
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D.  SURFACE  WAVES 

If  <u  tends  to  infinity  as  k  does,  we  must  proceed  differently. 

2 

Then  v  - 1  and  (5.17)  becomes 

tanh  ks  [l-tanh(ky1+o)tanh  k(y2+h>]  ~  tanhCky^+cO-tanh  k(y2+h)  (5.36) 
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Thus 


1  -  ranh  (ky^+of)  ~  tanh  (ky^+a)  -  1  (5.37) 

It  follows  that  a  must  tend  to  +»  as  k  does.  Then  (5.10)  yields 

u)  ~  (kg)^  (5.38) 

This  is  the  result  for  the  surface  wave. 

E.  DEEP  OCEAN'S 

In  the  oceans  the  depth  is  so  large  that  kh  »  1  even  for  the 
smallest  practical  value  of  k.  In  this  case  it  is  possible  to  simplify 
some  of  the  preceding  results.  For  example,  in  (5.17),  we  can  set 
tanh  k(y2+h)  =  1  +  ...  .  Then  the  solution  of  (5.17)  for  small  k  can 
be  carried  beyond  the  leading  term  given  in  (5.29)  and  (5.30)  with  the 
result 

u>  =  (g/|y1|  )^k[i  -  Xpl-  k  +  . ..]  ,  n  =  0  (5.39) 

(  '  ■& 

uj  =  -iSJLL  k[l  -  -§-*  k  +  ...]  ,  n  »  1  (5.40) 

n  tT 

In  writing  these  results  we  have  assumed  that  cu^/g|y^|  «  1,  since  in 
the  oceans  this  number  is  typically  of  the  order  IQ-3. 


VI.  WAVE  HEIGHT  DUE  TO  A  MOVING  SOURCE  OR  DIPOLE 


The  wave  height  7)  (x,z,t)  has  been  determined  for  a  unit 

OUUi uL 

point  source  of  fluid  moving  with  the  constant  velocity  -v  at  the 


depth  yQ.  The  asymptotic  form  of  T) 


source 


far  from  the  source  is 


t]  (x+vt,z) 
source  5 
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Yj — 5 — TTSk  cos  Ckr-(UT+TT/4) 
(v  -c2)V 


(6.1) 


(Ref.  1,  Eq.  11.14).  The  sum  in  (1)  is  over  the  modes  of  wave  propaga¬ 
tion,  and  for  each  mode  over  the  roots  k  and  t  of  (2.7)  and  (2.8).  The 
functions  w(k)  and  c(k)  =  u>/k  are  determined  for  each  mode  as  in  Sec¬ 
tion  5  and  B  is  given  by  Ref.  1,  (10.9). 


,  2  A.-1  C»2-N2(0)]w(yo) 

d  =  -  I — «r — I  yj  — o  ~  o  — 1 "  * 

(m>  O  -N^(o)]wk(o)  -gwky  (o)] 


(6.2) 


Here  w(y)  =  pD(y)  v(y)  where  v(y)  is  a  nontrivial  solution  of  (5.1) 
and  (5.3). 

If  the  source  is  a  dipole  of  unit  strength  oriented  along  its 
direction  of  motion,  the  wave  height  can  be  obtained  by  differ¬ 

entiating  (6.1)  with  respect  to  -x.  Only  the  phase  -0  of  the  cosine 
need  be  differentiated  and  in  view  of  (2.5)  its  derivative  is  -kr  + 

X 

wt  .  Alternatively  we  can  obtain  -0  from  (2.12).  In  either  way  we 

X  X 

obtain 


Vpole(x+vt>z>  ~-E- 


k(c-c  )  B 


v(l-ccgv-2)(v2-c2)^ 


sin  (kr-uiT-t-^)  (6.3) 
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VII.  ILLUSTRATION 


A.  NORMALIZATION 

Figures  (1)  and  (2)  are  drawn  for  a  density  distribution  consist - 

2 

ing  of  a  "thermocline"  of  thickness  s  =  y,-y0  within  which  N  is  con- 

2  * 

stant,  and  outside  which  N  =  0  (i.e.,  the  case  treated  in  Section  V). 
We  collect  the  dimensionless  formulae  used  in  the  construction.  All 
distances  are  normalized  with  respect  to  the  thermocline  thickness  s: 

X  =  x/s;,  Y  =  y/s,  Z  =  z/s;  X  =  ks  (7.1) 

ii  n  n 

and  all  frequencies  to  the  Vaisala  frequency  N,  so  that 

0  =  u>/N,  C  =  O./K  =  c/(Ns),  V  =  v/ (Ns)  (7.2) 


F  =  v/c  =  VK/Q  is  a  "Froude  Number,”  measuring  source  velocity  rela¬ 
tive  to  the  speed  of  internal  waves.  We  consider  only  F  >  1. 

The  dimensionless  amplitude  along  any  line  of  constant  phase 
0  =  0,  2tt,  ...  is  given  by 


where 


T)  ..  L 

'source  2" 2  B 


Ns5"*72  Z^V 


B 


^dipole  _  2"^  B0(1-F“2)"^ 


_  ,2a,*  _  N2Q  j*0%My0)\ 

vn  ;  g  cosh  ot  )  pQ(o)  j 


(6.4J) 

(6.5*  > 

(6.6») 
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The  expression  P0(y0)v(y0)/P0(°)  in  (6.6!)  depends  on  the  depth  of  the 

source,  yQ  =  sYq,  in  accordance  with  (5.5)  -  (5.8).  The  simplest  case 

is  that  of  the  source  above  the  thermocline,  Y^  ^  Y.  in  which  case  we 

o  1 

find 

po(yo)v(yo)  N2s  n2 

~-p-"(0)  '  '  =  sinh(KY0+&)’  tanh  01  =  ”cT  IT  (7’3) 


B.  NEAR  WAKE 


Equation  (4.6)  with  determined  by  (5.33)  can  be  written  in  the 
dimensionless  form 


Z 


A 

nn 


3 

2(X-0V)  |  2 
30V  1 


(7.4) 


Here  0  =  2tt,  4rr,...,  and  n  =  1,2  ...  designate  successive  crests  for 
various  modes.  The  cases  0=0  and  n  =  0  are  beyond  the  scope  of  the 
present  approximation.  The  wave  crest  can  now  be  constructed  for  any- 
specified  n  and  0. 

For  given  n,  0,  X,  and  Z,  the  amplitudes  can  be  obtained  as  follows. 
First  we  eliminate  between  (4.3)  and  (4.5)  to  obtain  z  =  cg0/N  +  ...  , 

with  0^  determined  by  (5.35).  This  leads  to  the  dimensionless  formula 


Z  ~  0 ( nn ) 2  (K+2)3  (7.5) 

from  which  K  can  be  calculated.  Furthermore  from  (5.33)  and  from  the 
definition  of  F  we  have 


Q  ~  1  - 


(7.6) 


F  =  VK/n 


(7.7) 
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By  differentiating  (7.6)  we  find 


~  3(mr)"4  (K+2)5 


(7.8) 


C.  FAR  WAKE 

Equation  (3.6)  can  now  be  written 


Z  =  (F^  -  l)"kx-2Xo^  +  ...] 


(7.9) 


where  F^  and  Xq  are  given  by 


F  = 


v 


o  a. 


(7.10) 


0ia2v' 


x° =  4(v2:»^ 


(7.11) 


The  mode  n  =  0  corresponds  essentially  to  a  thermocline  displacement, 
and  the  modes  n  =  1,2,...  to  thermocline  distortions.  We  need  to  treat 
these  cases  separately. 

From  (3.1),  (5.39)  and  (5.40)  we  find 

=  (g 7 1 y-j^i  )^5  u>2  =  (9,|y1l)^!y1l/2  for  n  =  0  (7.12) 

<u.  =  (g's^/nrr,  =  (g ,s)^s(nTT)-3  for  n  »  1  (7.13) 

By  using  (7.12)  and  (7.13)  in  (7.10)  we  obtain 

Fq  =  for  n  =  0;  Fq  =  nnV  for  n  »  1  (7.14) 

We  now  define  6(n)  by 


6(n)  =  1  for  n  »  1 
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6(n)  =  2  for  n  =  0; 


(7.15) 


i^fWWW'JM JXWMW  I'.IIU^H  l  \Kl<im  W*',  tM  WW'A'W-.  wmW4JW.!!«!llU  JJi'iJfi,KUU»M'  WW'JW  -  iJGgWPU  -  ■A^  W„U4*  i  M* 


^  >WO\-'-,-  3FJ 


l\tsctv; 


-^.  '•-  r-..  K[..r-r**7fl.“\*j» 


Then  from  (7.11)  -  (7.15)  we  obtain 

0V2F 

X  ^=  - £—  (7.16) 

0  6(n)(F^-l) 


We  now  use  (7.12)  and  (7.13)  in  (3.5)  to  get 

60F  (F2-l)^ 
<7-00 

^  2“i2 

vYT 


(7.17) 


From  (7.17)  we  can  find  K  in  terms  of  Z  and  0.  Then  from  (3.39)  and 
(5.40)  we  find 


Q  =  (V/Fq)K  -  6_1(V/Fo)3K2  +  ...  (7.18) 

Differentiation  of  (7. 18)  yields 

Km  ~  26~1£1  ~  26"1(V/Fo)2K]"3  (7.13) 


D.  RESULTS 

We  have  computed  wakes  for  the  following  cases: 
Source  depth  yQ  =  30  m 


Thermocline  depth 
Thermocline  thickness 

tt  tt  !t 

Vaisala  frequency 


y^  =  50  m 
s  =  10  m 
N  =  10” 2  sec- 


Figure  (1)  and  (2)  portray  the  near  wake  for  the  cases  V  =  VIS  and  10, 
corresponding  to  source  velocities  v  =  VNs  =  0.316  m/s  and  1  m/s, 
respectively.  The  X-axis  extends  from  X  =  0  to  200,  corresponding  to 
2  km  full  scale;  the  horizontal  Z-axis  is  drawn  to  the  same  scale. 
With  increasing  v,  the  wake  field  is  rapidly  concentrated  along  the 
source  axis,  particularly  for  large  n  and  0.  We  have  (improperly) 
used  the  n  »  .i  approximation  for  the  cases  n  =  1,  2.  The  case  n  =  0 
is  beyond  the  scope  of  the  present  treatment  for  the  near  wake*  and 
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the  far  wake  is  off-scale  in  the  example  shown.  The  computed  source 
functions  diminish  rapidly  with  distance  from  the  source  axis.  Unlike 
the  case  of  a  surface  (Kelvin)  wake,  internal  sources  moving  at  quite 
moderate  velocities  through  typically  stratified  fluids  produce  in¬ 
ternal  wakes  that  are  sharply  concentrated  along  the  source  axis. 
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or  a  source  velocity 
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ABSTRACT 


The  internal  waves  produced  by  either  a  moving  body  or  the 
collapsing  wake  behind  a  moving  body  in  a  stratified  fluid  are  calcu¬ 
lated  asymptotically  (at  large  distances  behind  the  source)  on  the 
hypotheses  of  small  disturbances,  the  Boussinesq  approximation,  and 
the  slender-body  approximation  (the  transverse  dimensions  of  the  body 
and  wake  are  small  compared  with  the  wavelengths  of  the  significant 
internal  waves). 

Explicit  results  are  given  for  two,  complementary  models:  (a) 
a  constant-N  model,  in  which  the  density  gradient  is  constant  and  (b) 
a  thermocline  model,  in  which  the  density  gradient  peaks  sharply 

in  a  thin  layer  and  is  elsewhere  negligible.  The  internal-wave  spec¬ 
trum  is  continuous  in  (a)  and  discrete  in  (b);  however,  only  the  domi¬ 
nant  mode  is  included  in  the  explicit  results  given  for  (b). 

A  WKB  solution  also  is  given  for  a  thermocline  model.  This 
approximation  does  not  give  an  adequate  representation  of  the  dominant 
mode  but  does  provide  estimates  of  the  contributions  of  the  higher 
modes  that  are  neglected  in  the  thin-thermocline  model.  These  contri¬ 
butions  of  the  higher  modes  that  are  neglected  in  the  thin-thermocline 
model.  These  contributions  are  typically  negligible  relative  to  that 
of  the  dominant  mode  in  the  neighbourhood  of  the  maximum,  free-surface 
disturbance . 
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I.  INTRODUCTION 

We  consider  the  disturbance  generated  by  a  horizontally  moving 
source  in  an  incompressible,  inviscid,  vertically  stratified  fluid. 

This  disturbance  comprises  the  near  field,  which  dies  out  more  or  less 
rapidly  with  distance  from  the  source,  and  the  radiated  field  which 
consists  of  internal  gravity  waves.  We  focus  primarily  on  the  radiated 
field,  but  emphasize  that  there  may  be  situations  of  interest  in  which 
the  amplitude  of  the  near  field  is  not  small  compared  with  that  of  the 
radiated  field.  In  particular,  the  radiated  field  in  a  steady  flow 
(uniform  translation  of  the  source)  appears  only  in  the  lee  of  the 
source,  so  that  the  near  field  must  be  taken  into  account  in  calcu¬ 
lating  the  disturbance  forward  of,  or  directly  over,  the  source. 

The  appropriate  similarity  parameter  for  the  generation  of  inter¬ 
nal  waves  by  a  moving  source  is  the  reduced  frequency  (or  inverse 
Froude  number) 


ft  =  Nt/U  (l.D 

IT  tt  It 

where  N  is  a  characteristic  value  of  the  intrinsic  (or  Vaisala)  fre¬ 
quency  of  internal  waves  (see  Eq.  2.9  below),  t  is  a  characteristic 
length  of  the  source,  and  U  is  its  speed.  The  frequency  spectrum  of 
the  internal  waves  is  (0,  N).  The  intensity  is  typically  a  rapidly 

JIlaX 

increasing  function  of  ft  (and,  therefore,  a  decreasing  function  cf  U) 
for  ft  <  ft  ,  say,  where  ft  is  a  characteristic  value  of  ft,  of  order 
unity,  at  which  nonlinear  phenomena  intervene .  Internal-wave  genera¬ 
tion  is  weak  for  ft  »  ftc,  and  as  ft  -*  ®  (U  -*  0)  the  flow  tends  to  a 
plug  type,  in  which  a  horizontal  column  of  fluid  is  pushed  in  front 
of  the  body.  x 
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We  develop  the  equations  of  motion  in  Sec.  II  on  the  hypotheses 
of  small  disturbances  and  the  Boussinesq  approximation  (in  which  only 
the  buoyancy  effects  of  density  stratification  are  included,  the  iner¬ 
tial  effects  being  neglected).  We  obtain  formal  solutions  of  these 
equations  in  Sec.  Ill  with  the  aid  of  integral  transforms  and  special¬ 
ize  these  to  a  moving  dipole  (by  which  a  body  may  be  approximated  if 
Q  «  1)  in  Sec,  IV  and  to  a  slender,  collapsing  wake  (a  region  of 
stirred  fluid)  aft  of  a  moving  body  in  Sec.  V.  We  give  explicit  cal¬ 
culations  of  the  internal-wave  field  for  a  constant-N  model  in  Sec. 

VI  and  for  a  thermocline  model,  in  which  N  peaks  sharply  in  a  region 
of  limited  vertical  extent,  in  Sec.  VII  and  Sec.  VIII. 

The  constant-N  model  is  characterized  by  a  continuous  spectrum 
(since  we  assume  the  fluid  to  be  either  infinite  or  semiinfinite)  and 
may  be  representative  for  laboratory  configurations,  although  finite- 
depth  effects  could  be  important  in  such  configurations.  The  thermo¬ 
cline  model  is  characterized  by  a  discrete  spectrum  and  affords  a  more 
realistic  model  for  the  ocean;  we  give  explicit  results  only  for  the 
dominant  mode  on  the  hypothesis  that  the  thickness  of  the  thermocline 
is  small  compared  with  both  its  depth  and  the  wavelength.  We  give  a 
WKB  solution  for  the  thermocline  model  in  Sec.  IX.  This  solution  does 
not  given  an  adequate  description  of  the  dominant  mode  for  a  thin 
thermocline,  but  it  does  provide  adequate  estimates  for  the  higher 
modes. 

The  disturbance  produced  by  a  moving  body  has  been  calculated 
previously  by  Hudimac  (Ref.  1)  for  a  two-layer  model  of  the  ocean  and 
by  Xeller  and  Levy  (Ref.  2),  Lighthill  (unpublished  papers),  and  Mei 
(Ref.  6)  for  various  models.  There  is  a  close  analogy  between  two- 
dimensional,  time -dependent  disturbances  and  three-dimensional  dis¬ 
turbances  produced  by  a  uniformly  translating  source.  Keller  has 
obtained  results  similar  to  (but  more  general  and  less  explicit  than) 
those  reported  here.  Many  reports  from  Hydror.autics,  Inc.  also  deal 
with  the  problem,  both  experimentally  and  theoretically.  Nevertheless, 
it  appears  that  some  of  the  results  given  here  are  new.  Perhaps  the 
most  interesting  are  the  asymptotic  approximations  to  the  respective, 
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lateral  strains  produced  at  the  free  surface  by  the  displacement 
(dipole)  and  wake  (quadrupole)  effects  of  a  submarine  that  is  small 
compared  v/ith  the  length  of  the  internal  waves,  i.e.,  Q  «  1.  Thus 
we  have 


|Hdl  ~  0.4a2tb^(h+|d-h|)  4(Nh/U^xd“^)  (1.2) 

-  5  _  11 

\\\  -  0.8a4b  ?(h+|d-h|)  ^  (Nd/Nh)2(U/Kh)H<  (1.3) 


where  a,  d  and  U  are  the  radius,  length,  depth  and  speed  of  the 
submarine.-  d  and  h  are  the  thickness  and  depth  of  the  thermocline, 
and  are  the  intrinsic  frequencies  at  depths  d  and  h,  and  xd 
and  a  are  the  respective  distances  behind  the  submarine  and  the  plane 
in  which  its  wake  begins  to  collapse. 
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II.  EQUATIONS  OF  MOTION 


We  consider  small  disturbances  in  an  inviscid,  incompressible, 
Boussinesq  fluid  in  which  the  (hydrostatic)  equilibrium  distributions 
of  density  and  pressure  are  pQ(z)  and  pQ(z)  and  z  is  measured  positive 
upwards.  Invoking  the  requirement  that  particle  density  be  conserved 
and  linearizing  the  equations  of  motion,  we  obtain 

p  =  pQ(z-t)  *  P0(z)-p^(zH  (2.1) 

V-v  =  m,  (2.2) 

and 

Pom:  =  “  ^  '  9(0,  0,  pj,  (2.3) 

where  p  denotes  the  density,  t  the  vertical  displacement  of  a  particle, 
v  the  velocity,  m  the  source  strength  per  unit  volume,  and  p  the  pres¬ 
sure,  each  as  a  function  of  the  Cartesian  coordinates  (x,y,z)  and  the 
time  t;  letter  subscripts  denote  partial  differentiation,  and  the 
triplet  -}  denotes  the  Cartesian  components  of  a  vector.  We 

seek  a  solution  of  (2,l)-(2.3)  for  a  prescribed  source  density  that  is 
introduced  at  t  =  0,  an  initial  displacement  $  (x ,y,z),  and  an  initial 

velocity  v  (x,y,z). 

'N,0 

Let  0  be  a  potential  such  that 

p  =  pQ  -  po0t  (2.4a) 

and 

v  =  Vq  +  {0X>  0y >  *t3  (2.4b) 
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Substituting  (2.4b)  into  (2.2)  and  invoking  the  continuity  equation  for 


v  •  v 
~0 


=  o, 


we  obtain 


where 


A0  +  ^  t  =  m, 


(2.S) 


(2.6) 


(2.7) 


is  the  two-dimensional  Laplacian,  and  the  operators  and  Imply 
partial  differentiation  with  respect  to  x  and  y.  We  have  assumed  that 
m  =  0  at  t  =  0;  if  m  =  m  at  t  =  0,  we  need  only  replace  the  right- 
hand  side  of  (2.5)  by  mQ  and  m  by  m-mo  in  (2.6).  Substituting  (2.4a,b) 
into  (2.1)  and  the  z- component  of  (2.3)  [the  x-  and  y-components  of 
(2.3)  are  satisfied  identically  by  (2.4a,b)],  eliminating  0  through 
(2.6),  and  invoking  the  Boussinesq  approximation  (thereby  neglecting 
Pq  except  where  it  is  multipled  by  g),  we  obtain 


VZZtt 


+ 


r)A*  =  mzt, 


(2.8) 


where 


N2  =  N2(z)  =  -gPQ(z)/p0(z)  (2.9) 

t»  tt  tr 

is  the  square  of  the  intrinsic  (Vaisala)  frequency. 

We  seek  the  solution  of  (2.8)  for  the  initial  conditions  (which 
follow  from  our  definitions) 


♦  =  t0»  0  =  ft  =  m  =  0 

and  the  boundary  conditions 


(2.10) 
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0 1 »  !  ❖ !  <  00  ( M  ,  |y|  -  ®) 


(2.11a) 


and 


t  =  0  (z  =  0,  -D), 


(2.11b) 


corresponding  to  a  free  surface  at  z  =  0  (which  acts  approximately  as 
a  rigid  boundary  for  internal  waves)  and  a  rigid  bottom  at  z  -  -D. 


A  convenient  measure  of  the  disturbance  at  the  free  surface  is 
the  lateral  strain, 


t 


^yy-(x>y’0’T)dT> 


(2.12) 


which  plays  a  significant  role  in  calculating  the  interaction  between 
the  internal  waves  and  p^e-existing  surface  waves. 
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III. 


FORMAL  SOLUTION 


We  define  the  transforms 


*  =  »xV’  M  =  £Wm’  *o  =  We 


( 3 . la , b , c ) 


where 


£( 


CO  1“ 

)  =  J  e_at(  )dt,  X-1(  )  =  (2n ±)~1J  eat(  )do  (Bo  >  0) 


-1“ 


(3.2a,b) 


3x<  > 


=/ 


-1 

fx 


e"iQ,x(  )dx,  3V  (  )  +  (2n)V  eiaX(  )da,  (3.3a,b) 


and  similarly  for  3  ,  with  x  and  <y  replaced  by  y  and  8,  respectively. 
Transforming  (2.8)  and  invoking  (2.10)  and  (2.11a),  we  place  the  result 
in  the  form 


O2  -  X2)(*  -  a'1^)  =  o-1Mz  +  kW*  ,  (3.4) 

where 

X  =  K^L*(N/a)2  (RX  >  0)  (3.5) 

K  =  Jct2+&2  (K  £  0),  (3.6) 

and  implies  partial  differentiation  with  respect  to  z.  The  boundary 

Ct 

conditions  for  ♦  are  given  by  (2.11b). 

The  Green’s  function  for  (3.4)  and  (2.11b)  is  determined  by 
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(3.7) 


-  X2)G(Z,C)  =  6(z-C) 

and 

G(0,C)  =  G(-D,C)  =  0  (3.8) 

and  yields  the  formal  solution  (after  integrating  the  term  in  M  by 
parts) : 

Kz)-a-1t0(z)  =  -a'-j/Gc(z,OM(C)dC  +  K  V^/G(z,e)N2(C)*0(e)dC  (3.9) 

We  have  suppressed  the  explicit  dependence  of  the  transforms  on  a,  8 
and  ct;  the  integrals  are  over  the  domains  o.f  M  and  which  we  assume 
to  be  of  finite  extern:,  and  6  is  Dirac’s  delta  function. 

Transforming  (2.6),  we  obtain 

§(z)  =  K"2(atz.^0z  -  M),  (3.10) 

which  completes  the  reduction  of  the  formal  solution  to  the  determina¬ 
tion  of  the  Green’s  function  and  the  evaluation  of  inverse  transforms. 
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IV.  MOVING  DIPOLE 


We  now  consider  the  disturbance  produced  by  the  dipole 

m  =  UD3  6(x+Ut)6(y)5(z+d)  (t  >  0),  (4.1) 

A 

which  is  introduced  atx=y=t=0  and  z  =  -d  and  moves  along  the 
negative-x  axis  with  the  uniform  speed  U.  The  parameter  D  is  the 
dipole  moment  and  has  the  dimensions  of  volume  (see  below).  Trans¬ 
forming  (4.1)  in  accordance  with  (3.1b),  we  obtain 

M  =  UDia(a-iaU)-16(z+d)  (4.2) 


Substituting  (4.2)  into  (3.9)  and  assuming  the  fluid  to  be  initially 
undisturbed  (ijr  =  0),  we  obtain 

^  =  -UDiora  1(a-iofLl)  ^.(ZjOj^  _  _d  (4.3) 


The  asymptotic  limit  of  $  as  t  -  00  is  determined  by  the  pole  of 
idplace  transform  at  a  =  iarll  (c 
equations  of  motion),  which  yields 


"the  Laplace  transform  at  a  =  iaU  (corresponding  to  ~  Udx  in  the 


Z~\-  ~  DG1(z)elaUt  (t  -  »), 


(4.4) 


where 


G3 (“ ^  "  "GC(Z’^|c  =  -d,  a  =  -j 


la'll 


(4.5) 
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We  determine  the  behaviour  of  X,  qua  function  of  a,  in  (4.5)  from  the 
antecedent  requirement  that  RX  >  0  as  c  approaches  the  imaginary  axis 
from  the  right: 

X1  H  XU  =  ioU  =  k*A-W«>5  (|«l  >k)  (4.6a) 

=  Ua'Jtt  (\a\  <  k),  (4.6b) 

where  k  =  k(z)  =  N(z)/U  (4.7) 

(^  also  has  an  essential  singularity  at  a  =  0,  which  makes  no  contri¬ 
bution  to  the  wave  field,  and  branch  points  associated  with  the  branch 
points  of  X,  qua  function  cf  a,  which  contribute  transients  that  die 
out  at  least  as  rapidly  as  1/t . )  Taking  the  inverse  Fourier  transform 
of  (4.4),  we  obtain 

♦  ~  D3fx+Ut  ^y1  Gl^  5  td(x+Ut,y,z),  (4.8) 

where  the  subscript  d  implies  dipole . 

Substituting  (4.2)  and  (4.3)  into  (3.10)  and  proceeding  as  above, 
we  obtain  the  corresponding  result 

0  ~  UD3"JUt  3"1  {icy/IC2)[Gl2(z)-6(z+d)]} 

s  0d(x+Ut,y,z)  (4.9) 

Substituting  (4.9)  into  (2.12),  we  obtain 

11  ~  <9/k>2C1z) I ■  Hd(x+ut,y)  (4.10) 

We  apply  these  results  to:  (i)  small  bodies  of  characteristic 
length  a  and  arbitrary  shape  and  (ii)  slender  bodies  of  characteristic 
transverse  and  axial  lengths  a  and  l,  where,  by  hypothesis, 


39 


ka  h  N(-d)a/U  «  1 


C4.il) 


and 

e  =  a/'t  «  1  (4.12) 

We  add  that  a  slender  body  for  which  k£  «  1  is  also  small. 

The  solution  to  the  problem  of  a  small  body  moving  with  uniform 
speed  U  follows  from  the  fact  that  the  flow  in  the  neighbourhood  of 
the  body  is  locally  potential  (ka  «  1  implies  that  the  effects  of 
stratification  are  negligible  over  a  region  of  scale  a).  Invoking  the 
well-known  result  that  the  potential  flow  past  a  body  is  equivalent 
to  that  induced  by  a  dipole  at  distances  R  that  are  large  compared 
with  a,  we  may  match  the  potential-flow  solution  to  the  solution 

(4.8)  in  an  intermediate  region  a  «  R  «  1/k  and  then  use  (4.8)  and 

(4.9)  to  determine  the  far  field  (Rayleigh- scattering  approximation). 
The  dipole  moment  is  given  by  Lamb  (Sec.  121a,  Ref.  5). 

D  =  V+V,.  (ka  «  1),  (4.13) 

where  V  is  the  volume  of  the  body  and  P0V*  its  virtual  mass  with  re¬ 
spect  to  axial  translation  in  a  homogeneous  fluid  of  density  pQ(-d). 

The  solution  to  the  slender- body  problem  follows  by  analogy  with 
the  corresponding  problem  in  aerodynamics,  cft  Ward  (Ref.  7).  Omitting 
the  details,  we  obtain 


t (x,y,a,t)  ~yS(§){ijrd(x+Ut-S,y,c)/D}d5 

(ka  «  1,  a  «  l)  (4.14) 

and  an  analogous  result  for  0,  where  S(x)  is  the  cross-sectional  area 
of  the  body,  and  the  integral  extends  over  the  body.  We  remark  that 
(4.14)  reduces  to  if  k-t  «  1,  corresponding  to  the  fact  that  V*  « 

V  for  a  slender  body. 
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V.  COLLAPSE  OF  MIXING  REGION 


We  consider  next  the  collapse  of  a  small  [in  the  sense  of  (4.11)] 
mass  of  fluid  that  has  been  stirred — for  example,  by  turbulence— in 
such  a  way  as  to  conserve  its  mass  but  alter  its  potential  energy  with 
respect  to  a  horizontal  plane  through  its  original  center  of  gravity, 


say  z  =  -d.  Our  definition  of  d  then  implies 

Iff  (z+d)pQ(z)dV  =  0,  (5.1) 

conservation  of  mass  implies 

///Cpo(2-V-,’o(2)3dV  *  -/ ff  poVV  =  °>  (5-2) 

and  the  potential  energy  is  given  by 

■  •///  (z+d)p0(z-’)0)dV  4=  -g JJJ i z+d)p^(zModV  (5.3a) 

=  Qp  (-d)N2(-d),  (5.3b) 

o 

where  Q  is  the  quadrupole  moment  of  the  region. 


Considering  now  the  second  integral  in  (3.S),  we  expand  G(z,£) 
about  Q  =  -d  to  obtain 

yG(z,C)N2(C)*0(C)d£  =  G(z,  -d)yN2(C)^0(C)d£ 

+  Gc(z,  -d)y(C+d)N2(C)^o(C)dC  (5.4) 


and  reduce  (3.1c)  to 


^0(C)  * JJ $Q(x, y,C)dxdy  (5.5) 

by  virtue  of  our  assumption  that  the  dimensions  of  the  mass  are  small . 

Substituting  (5.5)  into  (5.4),  we  find  that  the  first  integral  on  the 

o 

right-hand  side  vanishes  while  the  second  reduces  to  ON  (_d)  by  virtue 
of  (5.2),  (5.3),  and  the  Boussinesq  approximation.  Substituting  the 
resulting  approximation  into  (3.9),  we  obtain 

S'-cT1^  =  QN2(-d)K2CT"3G^(z,  -d)  (5.6) 

We  apply  this  last  result  to  a  collapsing  wake  in  the  lee  of  a 
small  moving  obstacle  on  the  hypothesis  that  the  fluid  in  the  wake  is 
mixed,  and  perhaps  also  augmented  by  turbulent  entrainment,  over  a  dis¬ 
tance  xQ  behind  the  obstacle,  at  which  point  the  turbulent  wake  begins 
to  collapse  and  releases  the  potential  energy  UE'(xq)  per  unit  time. 

The  resulting,  asymptotic  (as  t  -  00 )  disturbance  then  is  given  by 

t 

WQ  ~  -Q'(xo)N2(-d )f  dT£"^rJUT_x  J^V^tx.-d)}  (t  -  »)  (5.7) 

Carrying  out  the  integration  with  respect  to  t  and  invoking  the  fact 
that  (as  in  Sec.  IV  above)  the  inverse-Laplace  transform  of  the  result 
is  dominated  by  the  pole  at  a  =  iorll,  we  obtain 

♦-♦o  ~  3y1t<>‘2/i‘>3>G1U» 

o 

=  t  (x+Ut-xQ,  y,z),  (5.8) 

where:  k  is  given  by  (4.7);  G1  is  given  by  (4.5);  Q/(xq)  is  the  cross- 
sectional  quadrupole  moment  of  the  wake,  is  defined  as  in  (5.3b),  and 
has  the  dimensions  of  (length)4;  the  subscript  q  implies  quadrupole . 
Similarly, 
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VI.  CONSTANT- N  MODEL 


We  new  consider  the  specific  model  of  a  fluid  in  which  N  (and, 
hence,  also  k)  is  constant.  This  is  a  realistic  model  for  those  lab¬ 
oratory  configurations  in  which  the  effects  of  lateral  boundaries  may 
be  neglected.  It  is  not  a  realistic  model  for  typical  oceanic  config¬ 
urations,  but  it  does  provide  an  extreme  complement  to  the  thermocline 
model  of  the  following  sections.  We  give  special  consideration  to  the 
limiting  case  D  =  °°,  which  is  appropriate  for  oceanic  applications. 

The  solution  of  (3.7)  and  (3.8)  is  given  by 


r-%  _  sinh(xz)sinh[x(C+D)] 

- TqiwmTr^ -  <-z  >  g) 


(6.1) 


wherein  z  and  G  must  be  interchanged  if  z  <  G-  We  observe  that  G  is 

2 

a  meromorphic  function  of  X  ,  and  therefore  of  each  of  or,  8  and  o, 
for  finite  D,  and  has  the  Fourier- series  representation 


G(z,G )  =  -2D£ 
n=l 


sin(nnz/D)sin(nTrG/D) 
(XD)2  +  (nn)2 


(6.2) 


We  consider  first  the  limiting  case  D  -  for  which  (6.1)  re¬ 
duces  to 

G(z,G)  =  X-1e^sinhXz  (z  >  G?  D  =  ®),  (6.3) 


which  has  the  branch  points  of  X,  qua  function  of  each  of  a,  8  and  a. 

Substituting  (6.3),  together  with  the  complementary  result  for  z  <  Q, 

into  (4.5),  we  place  the  result  in  the  fom 

-X, |z+d|  X.(z-d) 

G1(z)  =  ^e  sgn(z+d)  ,  (6.4) 
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where  is  given  by  (4,6).  We  may  interpret  (6.4)  in  terms  of  a 
source  at  z  =  -d  and  an  image  at  z  =  d.*  We  carry  out  a  detailed  anal¬ 
ysis  only  for  the  first  term  and  move  the  origin  to  z  =  -d  with  the 
implicit  understanding  that  z  must  be  replaced  by  z+d  and  the  image 
solution  incorporated  in  the  final  results.  In  brief,  we  consider  a 
(dipole  or  quadrupole)  source  at  the  origin  of  an  unbounded  fluid  in 
which  N  is  constant, 


G(z,C) 


(6.5) 


Similar  results  may  be  obtained  for  0^  and  ^  by  substituting  (6.6) 
into  (4.8),  (5.9)  and  (5.8),  respectively.  We  recall  that  x  now  is 
measured  in  a  reference  f ramie  moving  with  the  source  (x  replaces  x  + 

Ut  in  the  development  of  Secs.  II-IV  above)  and  that  0^  is  an  asympto¬ 
tic  solution  that  is  strictly  valid  only  for  kx  -  ®  (although  experi¬ 
ence  suggests  that  the  asymptotic  approximation  is  likely  to  be  quali¬ 
tatively  valid  for  only  moderately  large  values  of  kx,  say  kx  >  1). 


We  obtain  stationary-phase  approximations  to  0d,  fd,  0^  and 
in  the  appendix  to  this  analysis.  Introducing  the  spherical  polar 


coordinates  R,  9  and  9  according  to 

'-The  Tmage  term  in  (6.4)  also  may  be  expressed  as  +ke 


sgn(z-d) 


x  =  Rcos6,  r  =  (y^+2^)2  =  RsinQ,  y  =  rcoscp, 

z  =  rsincp  (0  <  6  <  rr,  0  <  y  <  2tt)  (6.9) 

and  letting  kR  -*  ra  with  6  and  cp  fixed,  we  find  that  xCo'jP)  has  two  (no) 
points  of  stationary  phase  if  9  <  (>)  reflecting  the  fact  that 
internal  gravity  waves  (for  which  the  group  velocity  exceeds  the  phase 
velocity)  appear  only  downstream  of  their  source  in  a  steady  flow. 
Substituting  the  resulting  approximations  into  (2.4b),  we  obtain  the 
velocity  fields 

v^l  ~  -(k  DU/2rrR)  ^cotGsin^ <.C0o^9+sin^9sin^cp)^sin(kRsinc?) 

(kR  -  *>,  0  <  8  <  *>n),  (6.10) 

3 

and  Vq  ~  (k^Q,U/2nR)  9csc^9  (cos^p+sin^sin^cp/^costkRsincp) 

(kR  -  0  <  e  <  Jjtt),  (6.11) 

where 

0,=  (-sin8,  cosGcoscp,  cosPsincp]  (6.12) 

is  the  unit  vector  in  the  direction  of  increasing  8;  both  v^  and 
are  asymptotically  transverse  to  a  spherical  surface  with  center  at 
R  =  0  (a  well  known  property  of  internal  gravity  waves). 

The  maximum  velocities  given  by  the  approximations  (6.10)  and 
(6.12)  are  achieved  in  the  neighbourhood  of  8  =  0;  however,  the  approx¬ 
imations  are  not  uniformly  valid  as  8  -*  0,  partially  in  consequence  of 
the  restriction  kr  »  1  (implicit  in  the  staionary-phase  approximation) 
and  partially  in  consequence  of  the  slender-body  approximation,  which 
does  not  give  an  adequate  description  of  the  interference  among  the 
shorter  waves  (which  are  especially  important  in  the  neighbourhood  of 
9=0)  that  originate  at  various  points  of  a  source  of  finite  cross 
section.  Assuming  r  «  x  in  (6.10)  and  (6.11),  but  imposing  the  re¬ 
striction  kr  »  1  (so  that  1/kR  «  9  «  1),  we  obtain 
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~  -(k2DU/2n) |y| zr~4r  sin(kxz/r)  (kx  »  kr  »  1),  (6.13) 

and 

v  ~  (k“’Q/U/2TT)x^jy|  ^r”^r  cos(kxz/r)  (kx  »  kr  »  1),  (6.14) 

whe  re 

r  =  {O,  y,  z] .  (6.15) 

The  corresponding  approximations  to  the  lateral  strains,  as  de¬ 
fined  by  (4.10)  and  (5.10),  are  (we  omit  the  details  but  emphasize 
that  the  results  calculated  from  0^  and  0^  have  been  doubled  to  incor¬ 
porate  the  effects  of  the  respective  image  solutions  at  the  free  sur¬ 
face) 


T)^  ~  (k2dD/rr)xj y |  ^r“^sin(kxd/r)  (kx  »  kr  »  1)  (6.16) 


and 


Tiq  ~  -(k^Q'/rOx^lyl  5r-9cos(kxd/r)  (kx  »  kr  »  1),  (6.17) 

p  o  J'  I  I 

wherein  r  =  (y  +  d  )2.  The  maxima  of  ^  and  T|  with  respect  to  (y| 

are  given  by 

T)^  =  (kD/87Td^)(kx)sin(kx/V2)  at  y  =  d  «  x  (6.18) 

and 

T)  =  -0.045(Q7TTd4)(kx)3cos(2kx/3)  at  y  =  &/5d  '<  x.  (6.19) 
The  loci  of  constant  phase  for  T).  and  Ti  are  hyperbolae,  corresponding 

Li  M 

to  the  intersections  of  the  conical,  stationary-phase  surfaces, 
kRsincp  =  x(«s)  -»s),  with  the  free  surface;  the  loci  corresponding  to 
the  approximations  of  (6.16)  and  (6.17)  are 

(kx/x)2  -  Cy/d)2  =  1.  (5.20) 
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It  does  not  appear  possible  to  obtain  a  simple,  asymptotic  approx¬ 
imation  (for  kx  >>  1)  to  (6.7)  that  is  uniformly  valid  with  respect  to 
kr;  however,  we  can  obtain  an  approximation  that  is  valid  at  y  =  0 
(although  still  suffering  from  the  aforementioned  deficiency  of  the 
slender-body  approximation)  by  first  evaluating  the  Fourier  integral 
over  8  [Erdelyi  et  al,  Ref.  8,  Sec  1.5(27)j,  whence 


0d  =  -(DU/ 4n  ) 


/ 


10fXv  r(  2,  2, 

e  KqLIcv  (y  +z  ) 


k2z2}^(k2-.2)^,  (6.21) 


where  the  real  part  of  the  radical  is  non-negative,  and  XQ  is  a  modi¬ 
fied  Bessel  function  of  the  second  kind.  Differentiating  (6.21)  twice 
with  respect  to  y,  integrating  with  respect  to  x,  setting  y  =  0  and 
z  =  d  (in  the  reference  frame  with  origin  at  the  source),  and  doubling 
the  result  to  incorporate  the  effect  of  the  image  solution,  we  obtain 


T|d(x,0)  =  (D/2n2d) j  | of | {d(a2-k2)^}elorxda.  (6.22) 

-CO 

The  dominant  contribution  to  the  integral  in  (6.22)  comes  from  the 
neighbourhood  of  a  =  k,  which  yields 


1)d  ~  -(D/d2)(2k/Ti'5x)^sin(kx-^TT)  (kx  »  1,  y  =  0). 


(6.23) 


Similarly,  we  obtain 

T]  ~  (6Q'/d4)(2kx/TT3)^sin(kx-^TT)  (kx  »  1,  y  -  0).  (6.24) 
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VII. 


THERMOCLINE  MODEL 


We  consider  nov;  the  discrete  spectrum  of  internal  waves  associ¬ 
ated  with  a  thermocline  model,  for  which  (by  definition) 

G  £  N2(z)  £  N2(-h)  =  (7.1) 


f  N2(z)dz  -  N2b  -T  Ms  =  g'  ,  (7.2) 

-D 

where  -h  is  the  vertical  coordinate  of  the  thennocline,  that  is  the 
plane  in  which  N(z)  achieves  its  maximum  value,  N^;  A p  is  the  total 
increase  in  density  across  the  thermocline  (Ap  «  p  by  hypothesis); 
and  g'  is  a  reduced  gravitational  acceleration.  Setting 

cr  =  iu»  (7.3) 


in  (3.5)  and  (3.7),  we  obtain 

U2+  (K/uj)2N2(z)  -  k2}g(z,C)  =  fi(z-C)  (7.4) 


2 

Invoking  the  assumptions,  (7.1)  and  (7.2)  above,  that  N  (z)  >  0 

2 

and  that  the  integral  of  N  (z)  is  bounded  (a  nontrivial  restriction  if 
D  =  “),  we  infer  from  Sturm-Liouville  theory  that  there  exists  a  dis¬ 
crete  set  of  eigenvalues,  say  K  ,  and  eigenfunctions,  say  f  (z),  that 
satisfy 


•fIt  would  be  more  conventional  to  regard  the  wave  speed,  c  =  w/K  , 
as  the  eigenvalue  for  the  Sturm-Liouville  problem,  but  werind 
it  more  convenient  for  the  subsequent  development  to  introduce  K 
as  the  eigenvalue  and  to  regard  both  <u  and  K  as  prescribed.  h 
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[^  +  (KnN/l»)S  -  <2!fn(z)  -  0,  (7.5) 

fn<0)  =  fn  ('D)  =  0>  0.6) 

and  0 

/< N/0,)2fmfndz  =  6mn  '  0.7) 

-D 

where  6mn  is  the  Kronecker  delta.  Expanding  G  in  the  f  in  the  usual 
way,  we  obtain 


G(z,C)  =H(K2-Kn2)"lf  (z)fn(C)  (7.8) 

n 

Substituting  (7.8)  into  (4.5),  we  obtain 

GjU)  =  -L(K2-Kn)"lfn(-d)fn(z)  (uJ  =  Ua)’  (7*9) 
where  f'(-d)  =  (df/dC)^_  d- 

Referring  to  Secs.  IV  and  V  above,  we  seek  the  far  field  (kx  »  1) 
of  a  moving  source.  Substituting  (7.9)  into  (4.8)-(4.10)  and  (5.8)- 
(5.10),  invoking  the  Fourier  integral 

3T1(s2  +  a2)-1  =  isa'V3^  (fca^O),  (7.10) 

and  setting  a>  ~  Ua}  v/e  obtain 

CO 

Ud>$d)  =  -(D/*tC/  y^c—n^lf^-d) 
n  -® 

•  l  (iall/K2)f^(z)[l-  j  Qf|_1Yne"^Q‘  ~Yn')  I  yi  ]  ,  fn(z)}da 

(7.11) 
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t«q,tq)  =  ik2(-d)Q'(xo)/4„!5: 


ana 


*{-Uf^(z),  (iK2/a)fn(z)}da  ( j y |  >0), 


(7.12) 


T1  =  -(D/4it)£ /  K"2f '(-d)f^(0)elaX(Yne“Ynlyl-|Qfje~lay')da5 
n  -® 

(7.13) 

CO 

Hq=  tk2(-d)Q'(x0)/4„)£  /(ivn/a3)f^(-cl){i;(0)eil'x-Vnly'dc, 

^  _  00 

(7.14) 


where 


Vn  =  C»2  -  k2)%  (*Vn  4  0) 


(7.15) 


and  f^(0)  =  dfn/dz  at  z  =  0. 
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VIII.  THIN  TBERMOCLINE  APPROXIMATION 


We  carry  the  development  of  the  preceding  section  further  for  a 

2 

thin  thermocli ne,  for  which  N  (z)  differs  significantly  from  zero  only 
in  a  small  neighbourhood  of  z  =  -h,  where  it  exhibits  a  single,  sharp 
peak.  We  also  neglect  bottom  effects  by  setting  D  =  This  model  is 
perhaps  more  realistic  than,  but  in  any  event  complements,  that  of  Sec. 
VI. 

The  dispersion  relation  for  the  dominant  mode  of  a  thin  thermo- 
cline  may  be  expressed  in  terms  of  the  thermocline  parameters  N^  and 
b,  as  defined  by  (7.1)  and  (7.2),  and  the  depth  of  the  thermocline  on 
the  basis  of  the  assumptions 


N^b/U  =  k^b  «  1  and  b/h  «  1.  (8.1a,b) 

Setting  N  =  0  for  jz+hj  »  b  and  invoking  the  boundary  conditions  (7.6) 
and  the  requirement  that  f(z)  be  continuous  across  z  =  -h  as  b  -  0,  we 
choose  the  solutions  above  and  below  the  theimocline  in  the  form 


-c  f-cschKh  sinhKz-i 
f(z)  =  fhl  eK(z+h)  ! 


(8.2) 


where  f,  =  f(-h).  Integrating  (7.5)  across  the  thermocline  and  remark- 
*  *  it  2  2 

ing  that  both  f  -K  f  and  N  vanish  except  in  the  immediate  neighbour¬ 
hood  of  z  =  -h,  where  f "  is  discontinuous,  f  =  f^,  and  the  integral  of 
N^  is  given  by  (7.2),  we  obtain 


-f  if. 


K2f+(K1N/w)2f}dz  =  {f# 


-h+e 

-h-e 


-h+e 

+  (K/u))"!^  J  N2dz  =  -K(cothKh  +  l)f 
—  h—  € 


h 


+  g'oy^V 


2 

Solving  (8.3)  for  u)  ,  we  obtain  the  dispersion  relation  [Ref.  9, 
(5.3.7)]* 

=  feg/K(l-e-2lCh)  (8.4a) 

-  g/hK2(i-Kh  +  ...)  (Kh  -  0)  (8.4b) 

~  ^g'K  (K'h  -  “).  (8.4c) 

Similarly 

j^l  -  |u)|/U  -  k^KbU-e'2*11)}^  (8.5a) 

-  kh(bh)^K(l  -  %Kh  +  ...)  (Kh  -  0)  (8.5b) 

~kh(^bK)^  (Kh  -  »)  (8.5c) 


and 


d  uj 


'gl  ~  d^ 


li  d|a-,l 
=  U  1 


-  v 


'  *"l-(l-2Kh)e*’2ICh 


dK 


(l-e-2Kh)^ 


) 

(8.6a) 


-  (g 'h)\l  -  K  ;  •»  ...)  (Kh  -  0)  (8.6b) 

~  ^(g  //2K)^  (Kh-®).  (8.6c) 


We  note  that  (8.6b)  and  (8.6c)  intersect  at  Kh  -  \  and  serve  as  rough 
approximations  (with  maximum  errors  of  20  percent)  for  Kh  <  \  and 


*The  preceding  derivation  is  an  abbreviated  form  of  a  technique  used 
by  Lighthxll,  (Ref.  10)  and  Drazin  &  Howard,  (Ref.  11).  This  tech¬ 
nique  also  may  be  applied  to  the  higher  modes,  but  the  results  are 
rather  unwieldy.  Moreover,  the  contribution  of  the  dominant  mode 
to  the  free-surface  disturbance  will  dominate  the  contributions  of 
the  higher  modes  if  (8.1a)  is  satisfied  (see  end  of  Sec.  IX  below). 
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Kk  >  respectively,  in  the  subsequent,  stationary-phase  approxima 
tions.  Substituting  (8.2)  into  (7.7),  we  obtain 

=  ^(l-e“2Kh).  (8.7) 

We  emphasize  that  (8. 4)- (8. 7)  hold  only  for  the  dominant  mode*  (n  =  1). 

We  use  the  approximations  of  (8.2)  and  (8.5)  to  obtain  asymptotic 
approximations  to  the  lateral,  free-surface  strains,  T)^  and  "H  ,  on  the 
basis  of  these  hypotheses:  (a)  the  contributions  of  the  higher  modes 
(n  ^  2)  are  negligible  compared  with  that  of  the  dominant  mode  (we 
omit  the  subscript  1  with  this  understanding)  and  (b)  |aj  «  K ,  for 
which  a  sufficient  condition  is  U  »  (g'h)^.  The  latter  hypothesis 
permits  the  approximation 


Y  4  iK  (0  s’  a  «  K)  (8.8) 

in  place  of  (7.15)  and  the  neglect  of  |<*|e~lay^  compared  with 
in  (7.13).  Invoking  these  approximations,  substituting  (8.2)  and 
(8.7)  into  (7.13),  setting  z  =  0,  and  choosing  k,  rather  than  a,  as 
the  variable  of  integration  (thereby  regarding  a  as  the  eigenvalue 
for  prescribed  K  in  the  Sturm-Liouville  problem),  we  obtain 

CO 

T)d  ~  -(D/2tt iK2(da/dK)D('K)e“KH+i(ax"Kly!  )d:<,  (8.9) 

0 


where 


and 


DOO  = 


H  = 


-2Khx-l.,  -2^dN 

(1-e  )  (1+e  ) 

-1 

h  +  ) d-h | . 


(d  <  h) 


(8.10) 

(8.11) 


*The  oscillations  of  fn(z)  across  the  thermocline  do  not  permit 
the  approximation  i(z)  4  f.  for  n  >  1  in  the  integrands  of  (7.7) 
and  (8.3)  n 
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We  obtain  the  corresponding  approximation  to  1)  by  replacing  DK  by 
k2(0)Q/(Xg)K4/iQ'J  in  (8.9).  The  discontinuity  at  d  =  h  corresponds 
to  the  discontinuity  in  f'(z )  at  z  =  -h  and  is  an  intrinsic  character¬ 
istic  of  the  thin-thermocline  approximation. 

Carrying  out  a  stationary-phase  approximation  to  (8.9),  we 

obtain 


T)  ~  -D(2rrx)-^{'<2(d»/dK)jd2a/dK2|-^D(K)e-KH+i(ax-Ki>'i+  I  (8.12) 


where  K  is  determined  by 

da/dK  =  Cg/U  =  | y |/x  (H  «  jy|  <  (gh'^T3*).  (8.13) 

There  is  no  point  of  stationary  phase,  and  T)  is  0(x_1)  rather  than 
0(x  ^)  as  x  -  00 ,  if  | y |  >  (gh')2x/U.  A  saddle-point,  rather  than  a 
stationary-phase,  approximation  must  be  used  if  |yj/H  is  not  large; 

K  then  must  be  determined  by  replacing  jyj  by  |y|  -  iH  in  (8.13)  and 
is  complex. 

The  maximum  value  of  |f)  |  corresponds  roughly  to  the  maximum  of 
2  u 

K  exp(-KH),  that  is  kh  4  2,  which  yields  a  value  of  K  that  increases 
from  1/h  to  2/h  as  d  increases  from  0  to  h  and  then  remains  at  2/d 
for  d  >  h.  We  may  refine  these  estimates,  at  least  for  d  <  4h, 
by  utilizing  the  asymptotic  approximation  (8.5c),  the  substitution  of 
which  into  (8.12)  on  the  assumption  that  K  is  real  (|yj  »  H)  yields 

!Hdl  ~  D(2T1Ux)^(^g,)^|D(K)ie“KH.  (8.14) 

We  find  that  the  maximum  value  of  (8.14)  occurs  at  Kh  =  1.0  for  d  «  h, 
Kh  =  2.1  for  d  =  h-,  and  Kh  =  2.25  for  d  >  h,  so  that  KH  =  2  provides 
an  adequate  basis  for  an  estimate,  namely  (we  take  D  r  l) 

!  \)  I  max  =  0-2Dg'*(Ux)-Vt  (d  <  4h)  (8.15) 

at 

lyj/x  -  (q'H)Vw  (d  <  4h).  (8.16a) 
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If  d  >  4h,  we  must  use  (8.6b)  in  place  of  (8.6c),  the  principal  effect 
of  which  is  to  replace  H~^4  by  d-2h-^  in  (8.15)  and  (8.16a)  by 

|y|/x  4  (g'h)^U-1[l-(2h/d)]  (d  >  4h).  (8.16b) 

The  counterparts  of  (8.12)  and  (8.14)  for  Tl  are 

T)  ~  k2(-d)Q/(x0)(2T7x)"^R{K4a“3(da/diC)jd2a/dtc2r^ 

.D(K)e“KH+i(ax“Klyl^TT)}  (8.17) 

and  |l)|  ~  [N(-d)/Nh]2Q/(x0)  (U/nx) V1(8/g/)V1/4D(K)e"lCH.  (8.18) 

The  maximum  value  of  |T|  |  occurs  at  *H  4  11/4,  where  the  deviation 
of  D(K)  from  unity  is  small,  so  that 

|  Tjq  |  max  =  1.0[N(-d)/Nh]2Q,(xf))g  (8.19) 

at  |y|/x  4  0.2(g/H)^/U  (d  <  5.5h)  (8.20a) 

or  jy S/x  -  (g /h)^U~1[l-2 . 75(h/d ) ]  (d  >  5.5h).  (8.20b) 


We  also  note  that 


We  use  this  last  result  to  compare  the  lateral  surface  strains 

produced  by  the  dipole  effect  of  a  small,  prolate  ellipsoid  of  radius 

a  and  length  £  and  its  wake  on  the  hypothesis  that  the  wake  is  (or 

has  the  same  potential  energy  as  a  wake  that  is)  fully  mixed  and  of 

2  /  4 

radius  a;  then  D  4  2na  1/3  and  Q,  =  na  /4.  Substituting  these  results 
into  (8.21),  we  obtain 
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(8.22a) 


(8.22b) 


The  factors  a^/b-t  might  lie  between  10“  ^  and  10  ~  for  a  typical  sub- 

2  2 

marine,  U/c  might  lie  between  10  and  10  ,  and  [N(-d)/N,  ]  is  less 
g  o  11 

than  unity  and  might  be  as  small  as  10  if  the  submarine  is  well  out¬ 
side  of  the  thermocline .  It  follows  that ,  within  the  limitations  of 
the  hypotheses  implicit  in  our  model,  the  dipole  effect  is  likely  to 
dominate  the  wake  effect.  Both  effects  achieve  their  maxima  if  the 
submarine  is  in  the  thermocline  (d  t  h  =*=  H)  and  fall  off  rapidly  with 
increasing  d/h. 
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IX.  WKB  APPROXIMATION 

The  WKB  approximation  to  the  solution  of  the  Sturm-Liouville  prob¬ 
lem  posed  by  (7.5)  and  (7.6)  may  be  expected  to  yield  qualitatively 

accurate  results  for  all  but  the  dominant  mode  (n  =  1  below),  although 

2 

the  implicit  assumption  that  N  (z)  is  a  slowly  varying  function  ren¬ 
ders  it  quantitatively  accurate  only  for  those  modes  for  which  K  b 
»  1.  It  is  generally  inadequate  for  even  a  qualitative  description 

of  the  dominant  mode  of  a  thin  thermocline ,  for  which  Kb  «  1  and 
2 

N  (z)  varies  rapidly  near  z  =  -h.  It  is  consistent  with  the  WKB 
approximation  to  neglect  the  effects  of  both  upper  and  lower  bounda¬ 
ries  (the  implicit  restrictions  are  *nh  »  1  and  | D-h|  »  1,  re¬ 
spectively;  the  violation  of  the  former  restriction  is  likely  to  be 
qualitatively  significant  only  for  the  dominant  mode,  while  the  latter 
restriction  is  almost  always  satisfied  in  a  real  ocean).  Bearing 
these  remarks  in  mind,  we  rewrite  the  Sturm-Liouville  problem  of  (7.5)- 
(7.7)  in  the  form 

+  K^w(z)}fn(z)  =  0;  (9.1) 

fn(±c°)  =  o,  (9.2) 

CO 

and  /  wfmfndz  =  6mn,  (9.3) 

_00 

where  w(z)  =  uT2N2(z)  -  1  (9.4) 

is  the  weighting  function,  and  is  the  eigenvalue.  The  results 
presented  in  (7.8)  through  (7,15)  remain  valid  for  this  revised 
formulation. 
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We  proceed  on  the  assumptions  that  (z)  satisfies  (7.1)  and  has 

only  a  single  peak  (N  =  at  z  =  -h)  and  that  |<u|  <  (waves  for 

which  jujj  >  are  not  propagated);  then  w(z)  has  only  two  zeros,  say 

z,  and  z  .  such  that 
*0  u 


w(z^)  =  w(zu)  =  0  (zt  <  -h  <  zu) 


(9.5a) 


and 


w(z)  >0  (z,  <  z  <  z  ) . 

f  -t  u 


(9.5b) 


We  also  define  arg  w2  =  0  for  w  >  0  and  infer  arg  w-2  =  1>tt  for  w  <  0 
from  the  requirement  fto  >  0  (or,  equivalently,  dtu  <  0)  and  the  facts 
that  N7(z^)  >  0  and  N/(zj  )  <0.  We  then  may  pose  the  WKB  phase  inte¬ 
gral  in  the  forms 


P(z)  =  j"  w^dz 

2, 


=  P^2U)  +  iQ^z)  (z  >  zu) 


>0  (z,  <  z  <  z  )  ,  (9.6) 


=  -  iQ^(z) 


(z  <  Z<t) 


where 


Q» ( 2 )  = 


■/ 


dz 


(9.7a) 


and 


z 

Q^z)  -J  (-w)^dz. 


u 


(9.7b) 


Invoking  the  fact  that  w  -  -1  outside  of  the  thermocline,  we  obtain 


0^(2),  QyCz)  ~  | z+hj  (jz+hj  »  b).  (9.8) 

The  WKB  solution  of  (9.1)  and  (9.2)  is  given  by  the  following  (we 
omit  the  details  but  note  that  the  problem  is  analogous  to  that  of  the 
harmonic  oscillator  in  quantum  mechanics) 
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We  proceed  on  the  hypothesis  that 

lal  «  |««nl  >  (9.14) 

by  virtue  of  which  we  may  approximate  (7.15)  by 

vn  *  i«n  ,  (9.1S) 

where  argvn  is  determined  by  the  requirements  Ryn  s  0  and  Aa  <  0,  and 
neglect  the  term  in  exp(-|ay|)  in  (7.11)  and  (7.13).  We  also  use  the 
approximations  of  (9.8)  outside  of  the  thermocline.  Invoking  these 
approximations,  substituting  (9.9)  into  (7.13),  and  restricting  the 
range  of  integration  to  that  of  the  propagated  waves  (|a|  <  k  ;  waves 
for  which  |a|  >  k  are  not  propagated  and  are  negligible  for  kx  »  1) 

T)  ~  -(D/2tt)£  A  ft/*k iK  2eiax-(2n-1)ilCi(lyl“iH)d0f  ’ 
a  n=l  n  ‘'o  1 

(h,  |d-h|  »  s)  ,  (9.16) 

where  An  =  (2n-l)[TT+(2n-l)"1]"1[sgn(h-d)]n  ,  (9.17) 

is  given  by  (9.13b),  and  H  is  given  by  (8.11).  Introducing  the 
change  of  variable 


a  =  k  sinC 


(9.18) 


and  the  parameter 


ixn  =  (2n-l)(kxs)"1(  |y|-iH)  (-Sjn  <  argp,  <  0)  , 


(9.19) 


we  rewrite  (9.16)  in  the  form 

JgTT 


T)  ~  -(kD/2TTs2)  Y^A  ft  f  isin2£sec^£e 
d  n=l  n  'l 


ikxsinC(l-nnsec  (9.20) 


* 


i 


'I 


3 


r 


*•  * 
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The  integrand  of  (S.20)  has  a  saddle  point  at  the  point  deter¬ 
mined  parametrically  by 

sinC  =  a/k  =  vn  nn  =  (l-v2)2(i+v2)-l  ,  (9.21a, b) 

where  is  a  complex  number.  The  contribution  of  this  point  dominates 
the  asymptotic  approximation  (as  kx  -*  ®)  to  the  integral  (after  an 
appropriate  deformation  of  the  path  of  integration)  if  |u  |  <1  and 
ll-M.^1  is  not  small.  Carrying  out  the  saddle-point  approximation,  we 
obtain 

\  ~  (kD/s2)^  Adn(kx,  **n^  ’ 

n=l 

3  _  3 

Adn  =  ^CTrkx)‘^Anft{v2  (1-v2)  2(l+v2)^(5+v2)“V(Xn"3TT/4)}  .  (9.22a,b) 

where  Xr  =  2kxv3(i+v2)-1  .  (9.23) 

Similarly,  starting  from  (7.14),  we  obtain 

OS 

Tlq  ~  [N(-d)/Nh32[Q,(x0)/s4]EiAqn(Hx,  nn)  , 

Aqn  =  J2(TTkx)'^(2n-l)2AnR{v^(l-v2)  l+^)^(3+v2)"^el(Xn"37T/4)} 

(9.24a,b) 

The  largest  terms  in  the  modal  summations  of  (9.22)  and  (9.24) 
are  those,  if  any,  for  which  | ix  |  is  small  and,  from  (9.21b), 

vn=  1  -  (^n)^+  0(Un|)  .  (9.25) 

Substituting  (9.25)  into  (9.23)  and  retaining  only  the  dominant  terms 
in  each  of  the  real  and  imaginary  parts,  we  obtain 

Xn  4=  kx  -  i(2n-l)^(kx/s)H(y2+H2)^jy|3^  =  kx  -  iXni  (9.26) 
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Substituting  (9.25)  into  (9.21)  and  (9.13b),  we  obtain 


ja/Kj  =  (2n-l)"\s|2nnl^  , 


(9.27) 


so  that  the  approximation  (9.25)  is  consistent  with  the  restriction 
(9.14). 


We  use  the  approximations  (9.25)  and  (9.26)  to  obtain  the  esti¬ 


mates 


| A,  1  =?  2  V^[n+(2n-l)"1]"1(2n-l)^(kx)V(yW)"8e"Ani  (9.28) 


U  V  O  O  —  Y 


and  I A  h2  T  TT"kTT+(2n-l)"-L]"-L(2n-l)‘f(kx)‘ts‘t(y^)  °e"Ani  . 

1  qn 1 

(9.29) 


v-l-i-1, 


5  5  7  7 

8 -X  . 


Assuming  |y|  »  H,  we  find  that  |Adn|  has  its  maximum  at 

1 y 1 /H  =  (2/9)(2n-l)(H/s)kx  »  1  (9-30) 

and  similarly  for  (A^J ,  with  2/9  replaced  by  2/49.  These  maxima  are 
fairly  sharp  (in  |y|/H)  and  therefore  can  be  achieved  by  only  a  single 
mode  at  any  given  point.  The  corresponding  maxima  in  |^dl  and  |^lql> 
neglecting  all  modes  except  that  for  which  (9.30)  and  its  counterpart 
for  |Adq|  are  satisfied,  are 

3 

h,|  =  0.025(2n~l)"^(kD/s^H2)(kx)"'5  (9.31) 

1  d'max 

7 

and  l^qlmax  =  0-15(2n-l)“^[N(-d)/Nh]2[Q/(x0)/s^H2] (kx)~^  (9.32) 

Comparing  these  maxima  with  those  of  (8.15)  and  (8.19)  for  b  =  4s/3 
and  n  =  2  (typically  the  most  important  of  the  higher  modes),  we 


SSBKS@3SeSEGSK5KS33B  &3T.V 
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conclude  that  the  contributions  of  the  higher  modes,  relative  to  those 
of  the  dominant  mode,  to  |TL|  and  |T)  |  are  not  likely  to  exceed 
0.05(H/s)^^  and  0.12(s/H)3/^,  respectively.  The  former  ratio  could 
be  larger  than  unity,  but  only  for  H/s  such  that  all  contributions  to 
| TL |  would  be  very  small;  the  latter  ratio  is  certainly  small — typi¬ 
cally  between  10-1  and  10” 
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APPENDIX  TO  ANALYSIS  B 
STATIONARY-PHASE  APPROXIMATIONS 


The  integral  in  (A6)  has  two  points  of  stationary  phase  at 
a  =  a  =  ±  kxz/rR  =  ±  kcosGsincp  (A7) 

O 

2  2 

and  6  =  80  =  -a  yx/r  =  =f  kcos  8csc6sincpcoscp  (A8) 

o  S 

if  x  >  0  and  no  points  of  stationary  phase  if  x  <  0;  R,  0  and  <?  are 
polar  coordinates,  defined  by  (6.9)  above.  We  also  obtain 


X(a  ,0  )  =  ±  kzR/r  =  ±  kRsincp 

o  S 

(A9) 

Xoa  =  T  R^/'kr?  , 

(A10) 

and 

xg0  =  ±  (r7R/kx2z)(y2R2+z2r 2)“^ 

> 

(All) 

where 

f(-a, 

the  upper  and  lower  signs  correspond  to  a 
0)  is  the  complex  conjugate  of  f(a,  0),  we 

£  0. 
find 

Assuming  that 

that  the  con- 

tributions  of  the  two  points  to  the  stationary-phase  approximation  to 
I  are  complex  conjugates,  with  the  end  result 

I  ~  (kxz/nr3R2)(y2R2+z2r2)^[f(ac,  0c)eikzR/r}  (x  >  0)  (A12a) 

=  (k/ifR)cot9sincp(cos2cfH-sin29  sin2cp  A(fC«s,  6s)elkRslnt|>)  (x  >  0). 

(A12b) 

Comparing  (6.7)  and  the  corresponding  representations  of 
and  $  to  (Al),  we  obtain 

=  ^D{-U(a2+B2)”^(k2-a2)^  ,  sgnz}  (A13) 

and 

f  {0q5  tq)  =  %Lk2Qfa['U(a2+32)^(k2-a2)^  ,  -(a2+02)sgnz} .  (A14) 
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Substituting  (A13)  and  (A14)  into  (A12b),  we  obtain 


~  (kD/2nR){-U,cotesin<p}(cos29  +  sin20sin2cp)^ 

•cos(kRsincp)  (kR  -  ®,  0  <  6  <  ^tt)  (A15) 

and  ~  (k2Q1/2TTR){-Ucsccp,  cot0}secftcsc20 

3 

•(cos  <fHsin20sin29>2sin(kRsincp)  (kR  -  <*>,  0  <  0  <  %n) 

(A16) 

Substituting  (A15)  and  (A16)  into  (2.4b),  we  obtain  (6.10)  and  (6.11). 
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SUPPORTING  ANALYSIS  C 

INTERACTION  OF  INTERNAL  WAVES  AND  GRAVITY  WAVES 

J.  Hartle 

University  of  California)  Santa  Barbara 


F.  Zachariasen 

California  Institute  of  Technology 


We  wish  to  study  the  effects  on  rather  small  surface  waves,  and 
specifically  on  their  height  and  slope,  of  the  existence  of  internal 
waves  in  the  region  of  ocean  through  which  the  surface  waves  are  mov¬ 
ing.  Surface  waves  damp  exponentially  in  depth  in  a  distance  compara¬ 
ble  to  their  wavelength,  so  for  small  surface  waves,  we  can  assume  that 
the  dimensions  of  the  internal  wave  are  very  large  compared  to  those 
of  the  surface  waves,  and  therefore  the  internal  wave  can  be  well  rep¬ 
resented  by  a  horizontal  depth  independent  current.  We  describe  this 
with  a  velocity  field  U(x,y,t),  which  is  a  function  only  of  time,  of 
the  horizontal  coordinates  x,  y,  and  with  no  vertical  component.  Fur¬ 
thermore,  we  may  expect  the  times  and  horizontal  distances  over  which 
U  varies  to  be  much  greater  than  those  over  which  the  surface  waves  of 
interest  vary. 

We  shall  ignore  viscous  and  other  dissipative  effects  for  the  sur¬ 
face  waves;  that  is,  we  shall  assume  that  damping  is  unimportant  over 
horizontal  distances  comparable  to  the  region  occupied  by  the  internal 
wave.  Typically,  we  shall  be  interested  in  dimensions  of  surface  waves 
which  dissipate  in  distances  considerably  longer  than  that.  On  the 
other  hand,  since  the  size  of  the  effects  we  are  interested  in  will  be 

characterized  by  the  parameter  U/c  ,  where  c  is  the  group  velocity  of 

y  y 

the  surface  wave,  we  are  also  most  interested  in  slow  (i.e.,  short 
wavelength)  gravity  waves.  Yet  these  are  also  the  waves  that  dissipate 
most  quickly.  We  must  therefore  strike  a  balance  between  the  two  re¬ 
quirements  . 

Finally,  we  shall  assume  incompressible  irrotational  flow  in  the 
region  of  ocean  occupied  by  the  surface  waves.  Irrotational  flow  is 
described  by  a  velocity  potential  0  which  satisfies 

?20  =  0  (1) 
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There  are  three  boundary  conditions  on  the  solution  to  the  equation: 

(1)  At  the  surface,  the  vertical  velocity  of  the  fluid  d0/5z  is  the 
same  as  the  time  derivative  of  the  height  of  the  surface  h(x,y,t),  so 
that 


dt  =  -§t  +  (V0  ‘  7)h  =  -Sz  (2) 

and  Bernoulli’s  equation  relates  h  to  the  derivatives  of  0: 

■§£  +  h  =  -gh  .  (3) 

(2)  The  second  boundary  condition  is  the  assumption  that  at  large  depths 
u  approaches  the  imposed  velocity  U 

u(x,y,a.t) - ►  U(x,y,t) 

Z  -*  -o° 

or  if  $  is  the  velocity  potential  for  U  so  thac 


0(x,y,z,t) - -  §(x,y,t) 

z  -* 

(3)  The  final  boundary  condition  is  the  initial  condition  that  for 
times  far  in  the  past  the  imposed  flow  vanishes  and  the  wave  approaches 
a  freely  propagating  wave  0Q, 


The  effect  of  the  imposed  flow  on  the  propagation  of  surface  waves 
is  expected  to  be  small  since  the  velocity  of  the  surface  current  is, 
in  general,  much  smaller  than  that  of  the  surface  wave  in  open  sea 
conditions.  The  hydrodynamic  equations  adumbrated  above  may  therefore 
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be  expanded  in  powers  of  U  and  only  the  linear  term  retained.  We 
begin  by  writing 


0  --  5  +  0q  +  01  (5) 

where  0q  is  the  velocity  potential  the  initial  wave  would  have  had  in 
the  absence  of  the  surface  current  and  0^  is  the  correction  due  to  its 
presence. 

The  equations  governing  0^  are  obtained  by  substituting  Eq.  (5) 
in  Eqs.  (1)  through  (4)  and  retaining  only  terms  linear  in  and  $. 
Noting  that  §  and  0q  separately  satisfy  Eqs.  (1)  through  (4),  we  have 

V20  =0  (6) 

with  the  boundary  conditions 

01(x,y,z,t) - -  0  (7) 

t  -* 

01(x,y,z,t)  ►  0  (8) 

Z 


and 


+  g  _±  +  2  ~  [(V§  +  701)  .  70o]  +  \  (V01  +  v$)  .  V(V0oy 
+  70q  •  7[*0  •  (7$  +  701)]  =  0  . 


As  a  consequence  of  the  small  amplitude  assumption,  terms  which  are 

quadratic  in  0  are  negligible  in  comparison  with  those  linear  in  0  . 

°  2  o 

Further,  since  0q  is  a  small  correction  to  0Q,  terms  like  010Q  may  be 

neglected  in  comparison  with  0^. 

The  boundary  condition  at  the  surface  then  becomes 
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,ww~>v 
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dt 


1  + 
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-  2  [v$  •  70  ] 

ot  0 


2  4^ 


'V 


(3) 


|  The  problem  is  now  to  solve  Eq.  (6)  with  the  boundary  conditions  (7) 
[  to  (9).  The  general  solution  of  Laplace's  equation  (Eq.  (6))  which 
l  '  satisfies  the  boundary  condition  (0)  may  be  written 


01(x,Z,t) 


-f 


d2k 


(2tt) 


doj 

“3  e 


i(k*x  -  mt)  kz 

^  p 


a1(Js,w) 


Since  0  ,  the  unperturbed  flow,  also  satisfies  Eqs.  (6)  and  (8),  it 
o 

has  a  similar  Fourier  decomposition  with  Fourier  transform  ao(k,u>). 
Here  and  in  the  following,  x  will  mean  a  two-dimensional  vector  in 


the  xy  plane.  Equation  (9)  then  becomes 


I 

'■ 

it 


5, 

$ 

;( 


? 


i 


; 

: 

k 

l 


a  (k,u,)  =  2iw~2  F(k,tu) 

gk  -  u) 

where  F(k,u>)  is  the  Fourier  transform  of  (U*  V0O)Z=O  given  in  terms 

of  a  and  the  Fourier  transform  of  11  by 
o 

2 

F(k,w)  =  i  f d  -S  d-r  U(k-£,  tu-v)‘  3  aoCa5v)  .  (10) 

J  (2tt y 

This  solves  the  problem  of  determining  the  perturbation  <t>1  to  a 
small  amplitude  surface  wave  0Q  caused  by  an  arbitrary  surface  cur¬ 
rent  U.  The  quantity  of  chief  interest,  however,  is  not  the  velocity 
potential  0  but  rather  the  height  n.  If  we  write  6h  for  the  change 
in  height  caused  by  the  surface  current,  then  we  have  from  Eq.  (3) 

6h  =  -  ?  [-si  +  a  •  7*0]  • 

Here  we  have  retained  only  terms  linear  in  U,  have  neglected  terms  in 

02,  0  0-,,  and  have  used  the  fact  that  there  is  no  vertical  displace- 
o  o  1 

ment  from  the  velocity  potential  i  in  accordance  with  assumptions  (1) 
through  (4).  For  o01/3t  we  have 
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F(k,«)  eiC&^  -  W:) 


...N  i(k‘X  -  o>t ) 
F(k,cu)  e  ~  ~ 


=  -  2 


•u~  •  70o +  2  /ttt?  e 


a  k  duj  gk  F(k,uQ  ei(k*x  -  cut) 
(2tt)^  gk  -  u/ 


Thus  for  6h  we  have 

6h(x,t)  =  2h_  F(k,»)  ei(fe’^  -  “,t) 

(2n)s  i»2  -  gk 


I' 


or  in  terms  of  the  Fourier  transforms  of  0Q  and  U  (cf.  Eq.  (10)) 


eh(x,t)  =  ll^s.  +  i  U(k.3,  .-V)  .  a 

y  J  (2tt)  cu  -  qk  J  (2n) 


£ 

(2tt) 


gk*>  (2n)' 


ao(3  •  v)  e 


i(k*x  -  cut) 

(V  rv 


(11) 


In  general,  the  nonlinear  effects  on  the  unperturbed  wave  will 
be  of  the  order  of  magnitude 


(V0)2/|| 


max 

2tt\ 


typically  a  number  like  1/50.  These  effects  are  therefore  large  com¬ 
pared  with  the  effect  of  the  surface  current  which  might  optimisti¬ 
cally  be  of  the  order  of  several  percent.  For  a  calculation  of  the 
total  wave,  these  nonlinear  effects  cannot  be  neglected.  If,  however, 
one  is,  as  here,  mainly  interested  in  the  change  in  the  wave  structure 
due  to  the  surface  current  as  calculated  from  Eq.  (11),  then  the  change 
arising  from  the  linear  part  of  the  wave  will  be  larger  by  the  factor 
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(hm3N//2TT\)  than  that  coming  from  the  nonlinear  correction.  To  calcu- 

JUaX  w 

late  6h  from  Eq.  (31),  we  can  therefore  replace  a  with  the  value 

0  k 

appropriate  for  a  plane  wave  with  wave  vector  k^,  frequency  u>  =  (g  kQ)  , 
and  amplitude  A: 


aoCk,u>)  =  (2tt)3  A  6(2)(k-k))  6(co-««o)  . 


One  finds,  then,  that 


ShCx.t)  =  =-1^2  +  i  f  A  »  U(k-k  k  <&•*** 

9  J  (2„)3  o’ 

(12) 


Writing 


+® 

U()s-k0,  «-»,)  -/  UQs-J^,  t')  dt'  , 

-.00 


the  id  integration  can  be  perfomed  by  evaluating 
+°° 

5-iu>(t-t' ) 

: — 5 -  =  I(t-t')  . 

u>  -  gk 


?<£ 
y  2n 


(13) 


In  order  that  0^  vanish  for  large  negative  times  the  poles  in 
the  denominator  of  the  integrand  in  Eq.  (13)  must  be  displaced  slightly 
into  the  lower-half-complex-u>-plane.  One  then  has 


I(t-t') 


0 

i  j-g-i,^cg  (t-t-/)  giv1^  (t-t^-j 


t  <  t' 
t  >  t' 


Making  these  substitutions  in  Eq.  (12)  and  displacing  the  k  integra¬ 
tion  by  an  amount  k^,  one  has 

t 


^(Xjtjt7)  dt 1 


A  e 


ik 

~o 


•  x-io 
~  o 
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7C 


6h(x,t)  = 


”i xw, - w-pv 


where 

J(x,t,tO  -j  ^jVsIfe^l  e+i'"o(t-t')  k. 


ik‘X 

U(k,t  )  e  ~  ~ 


-i^|k+ko|  (t-t')  i^gjk+kj  (t-t') 

_e  -  e  J  • 


Now,  the  wavelength  of  the  internal  wave  is  much  longer  than  that  of 
the  surface  wave  and  U(k,t)  will  be  sharply  peaked  about  k=0.  We  may 
therefore  expand  ^bout  k  0. 

I -  r -  k*^, 

\g|js+fe^l  =  \9  K  U  +  \  -y2  +  — )  • 


A  4; 

Writing  c^  =  \  k^  (g/kQ)^  for  the  group  velocity  of  the  surface  wave, 
we  have 


^,0  +  j 


-*  [m  +  c  .  k]  k  •  U(k,t ') 
.2  L  o  ~  ~o  7 


ik  •  [x  -  c  (t-t'). 3  2iw  (t-t')  ik  •  [x  +  c  (t-t')] 

|*  /v  rv  Q  av  /v 


Since  the  surface  current  varies  slowly  in  time  as  well  as  space 
[assumption  (4)],  the  second  term  will  be  small  compared  to  the  first. 
Therefore 


%  *  -  c^t-t')] 


J(x,t,t')  =  O0  -  i  j ^  k^.  U(Jc,t')  e 


=  [tu  -  i  c  .  V]  k  ♦  U(x  -  c  (t-t '),  t  ')  . 


If  we  denote  the  height  of  the  unperturbed  wave  by  hQ(^t)>  sc  that 


i  A  “o 

h0<X>t>  =  -  -g—  0 


then  we  may  write 


£Vi  ~r>  *  />  , 

£h  -  ZJZ  .  /  (u-  /  r  i  v 

- d.  L-  i  ko  - 


%  £  •  7]  k  •  U[x-c  (t-t'),  t']. 

*  ~0  ^  ^  /-^r-  '  7 


The  height  of  a  wave  can  change  at  a  given  point  and  at  a  given  time 
not  only  because  the  amplitude  of  the  wave  changes,  but  also  because 
the  phase  changes.  Since  it  is  the  change  in  amplitude  which  is  of 
chief  interest  as  far  as  the  identification  of  the  current  is  con¬ 
cerned,  it  is  important  to  separate  these  two  effects. 

The  general  wavy  can  be  written 


h(x,t)  =  A(x,t)  eiXC~,t) 


where  A  and  X  are  real.  For  small  6A,  6X  perturbations  away  from  un¬ 
perturbed  values  A  and  X  ,  we  have 

o  o 


6h  _  6A  .  ,y 

F"  "  TT  +  1  6X 

o  o 


Thus 


r,  ,6hv  s )  a  A 

-  Ke  (r)  = - - i  ^  *  *f  JSo*  Utx-c  <t-t#),  t1 


and  6X  =  Im  <$)  =  -  -J  U(*-cg(t-t '),  t')  dt' 


The  first  term  in  the  amplitude  enhancement  is  an  instantaneous  effect 
and  very  small.  The  second  term  is  a  time  integrated  effect  and  de¬ 
pends  on  the  gradient  of  the  flow. 
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In  addition  to  the  change  in  amplitude,  the  change  in  wavelength, 
frequency,  and  mean  square  slope  are  also  of  interest.  The  change  in 
the  wave  number  %  and  the  change  in  the  frequency  6u)  may  be  obtained 
from  the  change  in  phase  6X  through  the  relations 


6k  =  V(6X)  ,  6uj  =  -  d(6X)/5t  . 

2  2 
The  mean  square  slope  in  m  is  given  by  the  time  average  of  (7h)  . 

The  time  average  is  taken  over  a  period  of  the  surface  wave,  a  time 

which  is  short  compared  with  the  characteristic  variation  time  of  the 

surface  current.  In  this  case  the  time  average  may  be  expressed  in 

terms  of  the  complex  waves  by 

O 

m  =  \  '7h  *•  '7h‘"‘  . 


In  terms  of  the  change  in  amplitude  and  phase,  one  can  then  easily 
find  for  the  change  in  mean  square  slope 


6m 


m 


-( 


,  6A 

2  7T 

o 


k 

~o 


?6X 


) 


or,  since 


k  •  V  6X  =  k  •  6k  =  k  6k  , 
~0  ~0  ~  o 


(14) 


this  can  be  written  as 


6ni  _  „  6(kA) 

T*  “  2  FI 

m  o  o 

o 


Inserting  expressions  for  6X  and  6A/A  into  Eq.  (14),  one  has,  finally. 


.2  k  •  U(x,t) 

~2 - c~ 

mo  9 


v. 

-  3  - 


-.00 


c  (t-t'),  t']  dt' 

~g 
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In  order  to  investigate  the  magnitude  of  possible  enhancements 
of  amplitude  and  slope,  let  us  assume  that  the  surface  flow  U  has  the 
form  of  a  wave  propagating  with  a  phase  velocity  C  in  a  direction 
specified  by  a  unit  vector  n: 


U(x,t)  =  U(n  •  x  -Ct) 

/n/  '*»'  ^ 


In  this  case  we  have 


k  •  U(x,t) 
6A  _  ~o  ~  ~ 

£  2  c 


-  \  (k 

'■'•O 


7)/ 


dt'  k  •  ll[n 


'^-£g  t)+(n  •  Cg-OtO 


and  a  similar  expression  for  the  mean  square  slope  change.  Suppose 
now  we  follow  a  crest  in  the  internal  wave  which  for  simplicity  we 
assume  to  occur  when  the  phase  of  U  vanishes.  Then  x,  is  related  to 
t  by 


x  =  C  t  n 


and  we  have 


6A  _ 


U(0) 


2  c 


^  (k 
~o 


L 


U[(C  -  c  •  n)(t-t')]  dt'  , 


where  U '  denotes  the  derivative  of  U  with  respect  to  its  argument. 

Surface  waves  which  have  these  components  of  the  group  velocity 
in  the  direction  of  propagation  of  the  internal  wave  equal  to  the 
internal  waves’s  phase  velocity  may  experience  a  large  enhancement 
from  the  second  term.  For  these  waves 


c 

~g 


A 


•  n 


=  C 


and  the  enhancements  in  slope  and  amplitude  are 


k  .  11(0) 
~o  ~ 


3  (£  •  n)  (£ 

~0  <■>» 


U(0))  T 


and 

£  •  U(0)  .  A  » 

f  =  0  2  Cg  -  i  (i0  •  n)  (k0  •  U(0»  T 

where  T  is  the  time  that  the  surface  and  internal  waves  have  been  in¬ 
teracting.  If  T  is  long  enough  this  time-integrated  effect  will  be 
appreciable  and  will  dominate  the  instantaneous  first  term. 


As  a  particular  example,  let  us  take  for  U  a  sine  wave  with  wave 
number  K,  frequency  fi,  and  the  surface  current  in  the  direction  of 
propagation . 

U  =  i  U  sin  (K  ♦  x  -  n  t)  . 


If  we  denote  by  9  the  angle  between  K  and  then  6A/A  and  6m2/m2 
may  be  written 


6m 


U  cos  0 
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-  3  (TKUq)  cos  8 


sin  T  K(Cg  cos  0  -  C)‘ 

T  K(c  cos  9  -  C) 

9 


and 


(15) 


6A  Uo  C0S  6 
F*  "  2  c 


-  I;  (TKUq)  cos  o 


sin  T  K(c  cos  9  -  C) 

£ 


T  K(Cg  cos  0  -  C)~ 


] 


The  most  favorable  case  is  for  waves  traveling  in  the  same  direction 
as  the  internal  waves  with  cg  =  C.  For  these  waves,  the  dominant- 
effect  is 


~  -  h  TKU  =  .  3  TXU  . 

0 
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For  waves  which  travel  at  an  angle  with  respect  to  the  internal  wave 

but  have  their  component  of  c  in  the  internal  wave  iirection  equal 

y  2 

to  C,  there  is  the  same  enhancement  decreased  by  a  factor  cos'B. 

Waves  whose  velocity  component  in  the  direction  of  internal  wave 

propagation  is  greater  or  less  than  C  will  eventually  pull  ahead  or 

lag  behind  the  internal  wave.  This  is  indicated  mathematically  by  the 

decrease  of  the  bracketed  factor  in  Eq.  (IS)  for  large  T  if  c  •  cos  0 

y 

r  C. 

There  will  thus  be  a  strong  time-integrated  amplitude  and  slope 
enhancement  for  the  special  class  of  waves  which  ride  along  with  the 
internal  wave.  This  effect  is  proportional  to  the  gradient  of  the 
surface  current  and  to  the  time  of  interaction.  This  time,  in  turn, 
will  at  best  be  the  minimum  of  the  characteristic  times  of  decrease 
of  the  internal  wave  and  the  surface  wave  due  to  dissipative  effects. 
If  the  lesser  of  these  times  is  long  enough,  there  may  be  an  appreci¬ 
able  enhancement. 
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I.  INTRODUCTION 

We  shall  be  concerned  with  the  scattering  of  electromagnetic 
radiation  from  the  ocean  surface,  with  emphasis  on  frequencies  in  the 
microwave  region.  In  order  to  simplify  our  calculations  we  have  cho¬ 
sen  to  study  a  model  in  which  the  electromagnetic  wave  is  taken  to  be 
a  scalar  field.  We  believe  that  such  a  model  brings  out  all  the 
essential  physical  phenomena  involved  in  the  scattering  process,  ex¬ 
cept  possibly  in  scattering  from  the  sea  at  low  angles  of  elevation. 
When  the  scattering  occurs  near  Brewster's  angle,  the  scattering  of 
vertically  polarized  waves  should  be  strongly  suppressed,  a  phenomenon 
which  our  scalar  model  cannot  reproduce.  In  any  case  our  approxima¬ 
tion  scheme  for  studying  scalar  waves  breaks  down  at  small  angles  be¬ 
cause  of  the  phenomena  of  shadowing  and  multiple  reflection. 

The  actual  extension  of  the  theory  presented  here  to  the  full 
problem  of  vector  electromagnetic  waves  is  perfectly  straightforward. 
The  resulting  formulas  for  polarization,  etc.  will  be  given  in  a  fu¬ 
ture,  more  detailed  report.  In  the  following  discussion  we  also 
ignore  Doppler  effects  arising  from  the  fact  that  the  ocean  surface 
is  constantly  in  motion.  These  effects,  which  are  not  believed  to  be 
important  in  the  present  context,  will  also  be  treated  in  the  later 
report . 
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II.  SCATTERING  FROM  THE  OCEAN  SURFACE 


At  any  instant  the  air-water  interface  is  given  by  the  surface 
z  =  h(x,y),  the  origin  being  so  chosen  that  if  we  average  over  many 
instants,  <h(x,y)>  =  0.  We  assume  that  the  illuminating  radar  pulse 
is  of  such  a  short  duration  that  during  the  time  it  is  actually  inci¬ 
dent  on  the  sea  surface,  h(x,y)  changes  by  a  negligible  amount.  Since 
radar  pulses  are,  typically,  microseconds  long,  this  approximation 
should  be  excellent.  Our  problem,  then,  divides  into  two  parts:  (a) 
given  a  plane  electromagnetic  wave,  of  wave  vector  k,  incident  on  a 
fixed  sea  surface  z  =  h(x,y),  find  the  wave  scattered  in  any  direction 
and  (b)  find  the  average  power  scattered  in  any  direction  where  the 
average  is  taken  over  many  values  of  h(x,y)  corresponding  either  to 
many  different  times  or  many  illuminated  patches  on  the  sea  surface. 
Since  we  are  neglecting  instantaneous  motions  of  the  sea  surface,  we 
can  say  nothing  about  possible  doppler  shifting  of  the  frequency  of 
the  scattered  wave. 

For  frequencies  in  the  microwave  region,  the  index  of  refraction 
of  sea  water  is  well  represented  by 

2 

n  =  n  (1  +  ia/ion  e  ) 

O  0  0 

where  nQ  =  80,  a  =  3mhos/meter,  and  is  the  dielectric  constant  of 
vacuum.  The  quantity  cr/um^  ^  =  IC*Vv  describes  the  relative  impor¬ 
tance  of  conduction  and  displacement  currents  in  the  equation  of  mo¬ 
tion.  At  microwave  frequencies,  v  -  10“^  cps,  and  the  imaginary  part 
of  the  index  of  refraction  is  totally  negligible.  Therefore  we  may 
safely  think  of  our  problem  as  that  of  computing  the  scattering  of 
electromagnetic  waves  from  the  interface  between  two  purely  dielectric 
media  with  indices  of  refraction  1  and  80  respectively. 
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Boundary  value  problems  of  this  kind  can,  unfortunately,  be 
solved  exactly  only  if  the  boundary  surface  is  quite  simple:  plane, 
elliptic,  etc.  Therefore  we  are  forced  to  resort  to  approximate 
methods.  In  general,  there  are  two  sorts  of  boundary  which  admit 
simple  approximate  solutions.  In  the  first  case,  suppose  that  the 
radius  of  curvature  of  the  boundary  surface  is  everywhere  large  com¬ 
pared  to  the  wavelength  of  the  incident  radiation.  We  can  then  apply 
geometrical  optics  to  compute  the  intensity  of  scattered  radiation. 

In  the  second  case,  suppose  that  the  deviation,  A,  of  the  surface 
from  one  which  has  a  known  solution  is  everywhere  small  compared  to 
the  radar  wave  length  \.  Then,  by  a  simple  form  of  perturbation  the¬ 
ory,  the  scattering  can  be  computed  correct  to  order  a/X  (we  will 
shortly  show  how  this  is  done).  Therefore  if  the  sea  surface  h(x,y) 
can  be  written  as  h  =  hQ  +  h^,  where  h^  is  everywhere  small  compared 
to  X,  a-d  the  radius  of  curvature  of  hQ  is  everywhere  large  compared 
to  X,  we  can  combine  the  above  two  approximation  methods  to  get  a 
decent  solution.  Whether  or  not  this  can  be  done  clearly  depends  on 
the  detailed  nature  of  t-'ne  sea  surface. 

At  any  instant  the  sea  surface,  h(x),  can  be  written  as  a  Fourier 
integral  h(x)  =  J  dk  a  (k)  e  .We  can  obviously  make  the  decompo¬ 
sition  h  =  h  +  h.,  where 
o  1 


hQ  =/  dk  a(K)  elR‘x 
k<k 

c 

hl=f  dK  a(K)  eiR*X 

k>kc 

and  k  is  for  the  moment  arbitrary.  We  can  then  shov;  that  if  k  is 
c  c 

properly  chosen  the  surface  hQ  has  a  mean  radius  of  curvature  which 
is  large  compared  with  the  radar  wavelength,  while  the  mean  magnitude 
of  h-^  is  everywhere  small  compared  to  the  radar  wavelength.  Since  the 
sea  surface  is  a  random  process  we  can  really  talk  only  about  the 
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mean  valuec  of  the  coefficients  a(£).  All  such  information  is  con¬ 
tained  in  the  correlation  function 


p(x)  =  (h(x)h(O)) 

-/*  R(k)  eiS-5 

which  has  been  determined  experimentally  to  have  the  form 


C  =  -  x  10'3. 

TT 


The  cutoff  k^  corresponds  to  the  gravity  waves  moving  with  the  wind 
velocity  and  the  cutoff  k£  corresponds  to  very  short  (say,  1mm)  cap¬ 
illary  waves.  The  corresponding  functions  for  the  surfaces  h  and  h., 

o  j. 

are  then 


P0(x)  =  <ho(x)hQ(0)>  =  c/dR  ^  eiR‘X 

1c 


2 

Pl(x)  =  (h^(x)h^(O))  =  C  f dk  - 


iT 


ikx 


The  mean  square  height  H  of  the  surface  h^  is  then 

k„ 


H 


=  Pl(0)  =  2nC  J  dk/k3  =  nC  (k'2-^2)  =  2x10' 3kJ 


O 

and  the  mean  square  radius  of  curvature,  R',  of  the  surface  hQ  is 
given  by 
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We  should  like  to  satisfy  simultaneously  the  conditions  H/X  «  1  and 
X/R  «  1,  where  X  is  the  radar  wavelength.  Recording  to  the  above 
equations  H/X  =  0.007Xp/X,  and  X/R  =  0.29  X/Xc,  so  thar  if  we  choose 
X^  —  6X j  we  have  H/X  ^  x/R  —  0.05.  With  such  small  values  for  the 
expansion  parameters,  we  feel  safe  in  computing  the  scattering  from 
the  surface  hQ  by  geometrical  optics  and  in  computing  the  extra  effect 
of  the  surface  h^  by  perturbation  theory.  We  emphasize  that  our 
ability  to  make  both  expansion  parameters  small  simultaneously  is  a 
stroke  of  good  luck  depending  on  the  detailed  statistical  structure 
of  ocean  waves.  It  probably  is  not  possible  for  other  sorts  of  ran¬ 
dom  surface. 

We  now  have  to  show  how  this  approximate  calculation  is  carried 
out  in  detail.  In  order  to  demonstrate  the  ideas  involved  we  study 
the  scattering  of  a  plane  scalar  wave  from  a  surface  h(x,y)  which  can 
be  decomposed  into  two  surfaces  h  and  h^  in  the  manner  just  described. 
Once  we  have  solved  this  problem  it  is  not  hard  to  fold  in  the  com¬ 
plications  due  to  the  vector  nature  of  the  electromagnetic  field. 

In  a  medium  of  varying  dielectric  constant  n(x),  we  assume  the 
wave  function  ^  to  satisfy 

+  k^  n(x)]  t(x)  =  0,  k  =  w/c. 

This  means  that  if  n(x)  has  a  discontinuity  on  a  surface,  then  <|r  and 
its  normal  derivative  must  be  continuous  across  that  surface.  In  the 
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case  at  hand,  n(x)  takes  on  either  of  two  constant  values  1,  or  n, 

the  dielectric  constant  of  seawater,  jumping  from  one  value  to  the 

other  at  the  surface  z  -  h(x,y).  We  then  want  to  find  the  solution 

ik*x 

to  this  equation  when  a  unit  amplitude  plane  wave,  =  e  , 
k  =  k(cos9,  0,  -sin9)  is  incident  on  the  sea  surface  from  above.  (See 
Fig.  1). 

Corresponding  to  the  division  of  the  surface  z  =  h(x,y)  into  a 
part,  hQ(x,y),  with  small  curvature,  and  a  part  with  small  amplitude, 
hi(x,y),  we  can  write  n(x)  =  nQ(x)  +  n^(x).  nQ(x)  takes  on  the  values 
n  and  1  and  describes  the  air-sea  interface  z  =  hQ(x,y).  n^(x)  takes 
on  the  values  0,  ±(n-l)  and  is  nonzero  only  in  a  small  region  around 
the  surface  z  =  hQ(x,y)  as  is  described  in  Fig-  2. 

Let  us  suppose  that  the  solution  to  the  scattering  problem  for 

the  surface  z  =  h  (x,y)  is  known  and  let  it  be  called  ijf  .  Let  us  also 
o  ^  o 

define  6 ^  by  ijr  =  $  +W  where  is  the  desired  solution  for  the  surface 
z  =  (hQ  +  h^)  (x,y).  We  can  combine  the  two  equations 

[v2  +  k2  nQ(x)]  =  0 

[v2  +  k2  (nQ(x)+n1(x))]  (t0+6t)  =  0 

to  give 

[v2  +  k2  nQ(x)]  6^  =  -k2n1(x)if(x)  . 

This  equation,  in  turn,  can  be  put  into  integral  form  if  we  introduce 
the  Green’s  function,  Go(x,x')>  which  is  a  solution  of 

(V2  +  k2  n  (x))  G  (x,xO  =  o(x-x')  . 

X  O  v 

This  allows  us  to  write 

6 Kx)=  -k2/ dx  '  Gq(x,x  ')  n1(x/)  Kx  ')  • 


-  *sss«;s^*<3SSgCBBa»aB  3WT353I 


Since  rt^  is  nonzero  only  in  a  volume  which  goes  to  zero  as  h^  goes  to 
zero,  6^  is  of  order  h^.  To  this  order,  therefore,  we  can  replace  1r, 
within  the  integral,  by  ^ : 

§Kx)  =  -k2  J  dx '  Gq(x,x/)  n1(x/)tQ(x/)  . 

To  the  same  order,  we  can  actually  replace  the  volume  integral 
by  a  surface  integral  over  the  surface  z  =  hQ(x,y).  The  equations 
of  motion  satisfied  by  %  and  Gq  imply  that  at  this  boundary  surface, 
both  these  functions  and  their  normal  derivatives  are  continuous. 
Therefore  the  effect  on  the  integral  of  their  variation  over  the  small 
volume  in  which  n^  is  nonzero  is  higher  than  first  order  in  h., .  We 
can  therefore  write 

SKx)  =  -k2(n-l) / dS  Go(x,x/(S))H1(x,(S))  *o(x'(S)) 

where  the  surface  integral  is  taker,  over  z  =  hQ(x,y),  x'(S)  is  the 

three-dimensional  position  vector  of  the  element  of  surface,  and  h^ 

is  the  normal  distance  between  the  surface  z  -  h  and  z  =  h  +h_, 

o  o  1 

(taken  positive  or  negative  according  as  z  =  h^h^  lies  above  or  be¬ 
low  z  =  hQ).  Finally,  it  is  convenient  to  convert  this  into  an  inte¬ 
gral  over  the  plane  surface  z  =  0,  taking  p  =  (x,y)  to  be  the  position 
vector  in  that  surface 

6*(x)  =  -k2(n-l) J  dp  h1(p)  G0(x,x(p))  *0(x(p))  , 

where  x(p)  =  (p,  hQ(p))  and  h.^  is  exactly  the  quantity  earlier  called 

h, ,  the  vertical  distance  between  the  two  surfaces  z  =  h  and  z  = 

«*•  o 

ho+hr  The  geometry  of  the  transformation  is  best  explained  by  Fig. 

3.  Therefore,  if  we  know  £  and  G  on  the  surface  z  =  h  ,  we  can 

co  o 

calculate  5<jf,  correct  to  order  h^. 

According  to  our  assumption,  the  radius  of  curvature  of  h  is 

o 

everywhere  large  compared  to  the  wavelength  of  the  illuminating  radi¬ 
ation,  so  that  scattering  can  safely  be  computed  via  geometrical  optics 
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In  particular,  we  need  to  know  <r  on  the  surface  z  =  hQ.  Geometrical 
optics  means  that,  if  we  neglect  multiple  scattering  and  shadowing, 
the  field  at  a  point  on  the  surface  can  be  computed  by  replacing  the 
curved  surface  by  its  local  tangent  plane  and  imagining  the  given 
incident  wave  to  be  scattering  from  it.  It  is  easy  to  show  that  if 
the  incident  wave  is  e1  *x,  the  total  field  at  a  plane  boundary  be¬ 
tween  regions  with  dielectric  constants  n  and  1  is 


2cosa _  eik*x 

cosa  +  //n-sin^a 


where  cosa  =  £*n,  n  being  the  unit  normal  to  the  boundary 


We  also  need  to  know  Go(x,x')  with  x'  on  the  surface.  From  the 
equation  satisfied  by  Gq,  it  is  clear  that  G^Xjx')  represents  the 
total  field  generated  at  x'  by  placing  a  unit  source  at  x,  given  the 
boundary  specified  by  no(x).  If  x'  is  near  the  surface,  and  if  we  ne¬ 
glect  multiple  scattering  and  shadowing,  the  geometrical  optics  approx¬ 
imation  to  Gq  is  gotten  by  replacing  the  curved  surface  by  its  local 
tangent  plane.  The  solution  for  G^  in  the  presence  of  a  plane  boundary 
is  well  known.  If  we  set  x  =  £'R,  take  x  on  the  boundary,  and  let 
R-*“,  it  becomes 


G0(k'R,x)  = 


;-ik&'*x 

~?rrE 


2cosa ' 


cosar 


+  jr 


.  2  7 

sxn  a 


a  being  the  angle  between  k  and  the  local  surf  'a  normal. 

We  now  can  write  down  our  expression  for  the  field  6f(x)  when  x 
is  far  away  from  the  surface: 


R6t(k'R) 


k  (n- 
4tt 


d?5  hx(p  )T(a  (p )  )T(a 7  (p )  )eik(k^ ' )  *x (p  > 


where 


T(a)  =  2cosa/(cosa  +  -/n-sin^a)  . 


We  are  particularly  interested  in  the  field  scattered  back  along  the 

A  f  A 

direction  of  the  incident  beam,  in  which  case  k  =  -k,  and 
RSK-kR)  =  — j  do  hx(p)  T2(a(p))  e2i^*x(p)  . 

It  turns  out  to  be  convenient,  for  purposes  of  computing  the 
average  backscattered  power,  to  recast  the  expression  for  6f  in  a 
slightly  different  form.  First  of  all,  we  note  that  the  reflection 
coefficient  is  a  function  of  cosa  =  -n*it  where  n  is  the  local  surface 
normal.  In  turn,  n  is  a  simple  function  of  V— hQ(p),  so  that  we  can 
write  T  =  T(v— hoCp)).  If  we  introduce  the  Fourier  decomposition  of 
h^Cp),  h^OE),  we  then  have 

M(-kR)  =  4r2  f AX  ho(p)eiC2E+?;)-'-2ks:in0  ho(S)] 

Since  the  surface  h  (p)  is  one  which  satisfies  the  criteria  of  geo¬ 
metric  optics,  we  can  evaluate  the  integral  over  p  by  the  method  of 
staionary  phase.  This  means  that  the  only  important  contributions 
come  from  those  points  p  where 

((2k-t^)-p  -2ksin9  hQ(p))  -  0 
or 

h  (p)  =  (2k+£)/2ksin6 . 

P  o 


Therefore  we  have 


f  dp  exp[i((2k+?-)*p-2ksineh0)]. 


The  virtue  of  this  expression  is  that  the  arguments  of  T  no  longer 
depend  on  the  specific  surface,  so  that  the  averaging  process  is 
simplified. 
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To  compute  the  backscattered  power  we  need  <|^  +  6^  j 2 > 5  where 
^  is  the  backscattered  field  from  the  surface  h  ,  and  the  average  is 
over  the  various  possible  forms  of  the  sea  surface.  Since  different 
Fourier  coefficients  of  the  sea  surface  are  statistically  independent, 
and  since  6^  depends  on  h^  while  ^  does  not,  the  cross  terms  of  the 
form  #6^  vanish  upon  taking  the  average.  Therefore,  the  average 
backscattered  power  is  the  sum  of  two  terms,  (l^l2)  and  (|6^|2), 
which  we  shall  compute  separately. 

It  is  convenient  to  define  a  scattering  cross-section  in  order 
to  eliminate  the  distance  of  the  observation  point  from  the  sea  surface. 

r\ 

The  energy  density  at  any  point  is  just  |$j  .  If  a  finite  patch  of 
sea  surface,  of  area  A,  is  illuminated  and  we  observe  at  x  =  i 'R,  R 
very  large,  then  all  the  energy  at  x  is  flowing  in  the  direction  k/. 

If  the  antenna  subtends  a  solid  angle  AO,  the  total  received  power  is 
then  |<|r|  R  AQ  and  the  received  power  per  unit  illuminated  area  is 
U|2R2  Afi/A.  We  shall  define  the  quantity  a  =  |ijr|2E2/A,  so  that  an¬ 
tenna  power  is  a  A  AQ. 

i  1 2  2 

Let  us  first  compute  =  <| 6^1  >R  /A.  If  we  make  the  standard 
assumptions  about  the  Gaussian  nature  of  the  sea  surface,  and  make  the 
definitions 

<hQ(x)ho(0)>  =  Co(x) 

<K,a)h*(r)>  =  6  a-zf) 


we  find  that 


-  (k2(r.-l)/4n)2 


/ 


d Z  p,  (.1) 


rp4  f  2k+-t  ^ 

^  2ksinS ' 


dr  exp[i(2k+^) •r-(2ksin9)2(Co(0)-Co(r))] 
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,  .„  +.>*=+-  uhpn  we  write  R  we  mean  R  —  CR-,»Ry) 
,we  we  should  now  remember  that  when  we  wiiue  _  x  y 

S^e.0).  The  largest  contribution  to  the  integral  over  r  comes  an 
the  neighborhood  of  r=0,  where  we  write 

e  (r)  =  h  (1-ar  +  ...) 


so  that 


J dr  exp[i(2k+l).r-(2ksine)2CCoC0)-Co(r))]  “ 

j  dr  exp[i(2R+l).r-(2k3ine)2a  r2]  = 

J2L _  cxpC - C 2k+£ ) 2/4a( 2khsin9 ) 2 ]  =  f(2k+:t) 


(2khsin9)‘ 


We  note  that  in  the  limit 
a  is  rather  small: 


imit  a  -  0,  f(2t»l)  -  (2n)2  6  <2fc+l).  In  fact 


1 

h2a  =  -Cl/2)v2po(0)  =  (C/2)  y  dk/k2  =  "C  loge(kc/k1) 


=  2x10' 3  loga  0'wind/6,'radar) 

where  X  -  „  is  the  wavelength  of  those  ocean  waves  whose  velocity 
equals  S? wind  velocity.  For  a  10  m/s  wind  and  a  10  cm  radar  wave- 
length.  we  have  h2  .  0.9xl0'2.  This  means  that  in  terms  of  the  dime 
sionless  variable  |*M|  /  TksinO  the  width  of  f  is  about  0.2  in  a 
typical  situation.  This  width  decreases  slowly  with  wind  velocity. 

If  we  ignore  the  width  of  f,  replacing  it  with  a  delta  function, 

we  have  the  simple  formula 

0  =kVil!l4(0)p1(-2k). 

1  4  -*■ 
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If  we  include  the  effect  of  the  wind  broadening  of  f  we  see  that  a., 
is  proportional  to  the  average  value  of  T4p^(£)  ever  a  circle  of  radius 
=“  0.4ksin9  in  -t-space.  centered  about  X  =  ~2k  =  (-2kcos8,0).  In  the 
undistrubed  ocean  we  know  that  p-^Cw  =  2x10  l  /v.  and  therefore  could 
actually  perform  the  average,  if  necessary.  If  we  evaluate  the  zero- 
width  approximation  to  cr  in  the  limit  n-*°  (appropriate  for  the  sea 
surface  since  n  =  80),  we  get 

a1  =  (10"3/2n)  tan40 

for  the  backscattering  cross-section  from  the  undisturbed  sea. 

2  2 

We  now  muse  compute  a  =  }  )  R  /A.  Since  is  the  field 

generated  by  that  part  of  the  sea  surface  for  which  the  approximations 
of  geometrical  optics  are  correct  we  can  adopt  the  classical  results 
for  scattering  light  from  a  Gaussianly  rough  surface: 

=  (rrh2a  (2sin8 )4)-1exp(-eot29/ah2) 

for  backscattering.  We  note  that  ctq  falls  off  exponentially  as  8 

decreases  from  tt/2.  In  fact  we  can  easily  see  that  for  9  <  80°,  ct^ 

dominates  a  ,  while  for  8  ^  80°  the  reverse  is  true, 
o 

At  this  point  we  may  reasonably  summarize  our  results:  We  have 
found  two  basic  regimes  in  radar  backscattering;  one  occurs  when  the 
angle  of  elevation  is  large,  nearly  90°,  the  other  occurs  for  moder¬ 
ate  elevation  angles.  In  the  first  case,  the  backscattering  cross- 

2 

section  is  a  function  of  ah  ,  while  in  the  other  it  is  determined  by 
p^,  the  ocean  wave  power  spectrum,  at  some  appropriate  wave  number. 

The  quantity  ah2,  is  just  the  mean  square  slope  of  the  ocean  waves, 
which  in  turn  is  an  integral  over  the  complete  wave  power  spectrum. 
Therefore,  the  difference  between  the  two  regimes  is  that  in  one  case 
we  measure  p^,  at  a  point  in  wave  number  space  while  in  the  other  we 
measure  what  amounts  to  an  average  of  over  all  wave  number  space. 

This  distinction  will  turn  out  to  be  most  important  in  the  applications. 


Finally,  we  would  like  to  point  out  that  the  scattering  of  polar¬ 
ized  radiation  from  the  sea  surface  may  be  calculated  by  much  the  same 
methods,  although  the  formulas  are  much  more  complicated.  We  shall 
refrain  from  writing  them  down  here  since  nothing  essentially  new  in 
the  physics  of  radar  backscattering  is  introduced. 
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III.  THE  EFFECT  OF  INTERNAL  WAVES  ON  RADAR  BACKSCATTERING 


As  far  as  the  applications  considered  in  this  report  are  con¬ 
cerned,  we  need  to  know  the  effect  of  an  internal  wave,  over  and  above 
the  random  background,  on  the  radar  return  from  the  sea  surface.  The 
ocean  wave  heights  are  in  general  described  by  the  power  spectrum 
F({c)  where 


<h(x)h(0)>  = 

2x10” 

with  F(E)  =  - - -  for  the  standard  wind-generated  sea.  Hartle 

and  Zachariasen  have  shown  that  if  an  internal  wave  of  phase  velocity 
C,  wave  length  L,  and  maximum  surface  water  velocity  V  is  present, 
then  the  power  spectrum  is  changed  by 


yk  F(lc)  eiR'X 


6F(£)  „ 
F(£) 


T  2tt  V  cos^0 
o 

- E - 


sin(2Ti  TL”  (C  cos0-C) 
- __ - —2 - 

2tt  TL  (C  cos0-C) 

9 


where  0  is  the  angle  between  fc  and  the  direction  of  propagation  of  the 
internal  wave,  =  ^/g/k  is  the  group  velocity  of  the  surface  waves 
with  wave  number  ic,  and  T  ;  •  the  time  during  which  these  same  surface 
waves  have  been  acted  on  by  the  internal  wave.  It  is  convenient  to 
introduce  3"  =  TC/L,  which  is  just  T  measured  in  internal  wave  periods, 
and  e  =  Vq/C  so  that 


6f0c) 

F(k) 


2  sin(2nff(C  COS0/C-1)) 

e2,T:rcos  0  z^cok/c-n — 

y 


In  practical  cases,  e  turns  out  to  be  very  small.  Vie  note  that  for 
small  S',  6F/F  is  uniformly  distributed  over  k-space,  while  for  S'  »  1, 
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6F/F  is  very  sharply  concentrated  around  the  curve  Cgcos0  =  C.  This 
turns  out  to  have  a  profound  effect  on  the  size  of  the  radar  return 
for  different  values  of  3. 

We  recall  that  if  the  angle  of  elevation,  6,  of  the  radar  beam 
from  the  horizontal  is  not  too  large,  the  backscattered  power  is  pro¬ 
portional  to  the  average  of  F(R)  over  a  circle  in  R-space  centered  at 
R  =  (-2kradcos6,0)  and  with  radius  akrad_r,  where  a  depends  on  wind 
velocity  but  might  typically  be  0.1.  Let  us  suppose  that  we  have 
cleverly  chosen  k  ,  and  6  so  that  R  lies  on  the  curve  C  cos0  =  C. 
We  now  want  to  compute  the  ratio  6P/P  where  P  is  the  radar  return  from 
the  undisturbed  ocean  and  P+6P  is  the  radar  return  from  the  ocean  in 
the  presence  an  internal  wave.  Let  0  (|R-R0i)  be  1  for  |R-RQ|  < 

ak  ,  and  zero  otherwise.  Then 
radar 


6P  _ 
P~  " 


fak  6  F(R)  Q  (|k-k0|) 
jfdR  F(R)  Q  (|R-R0|) 


There  are  two  interesting  regimes  in  which  we  want  to  calculate 

this  ratio.  First  of  all,  if  2mT(C  cos0/C-l)  is  small  through  the 

9  2 

region  where  ft  is  non- zero,  we  have  6F  =  -e  2tt3cos  0F  and 


—  _e2n3'  cos^0  . 

o 

On  the  other  hand,  if  3"  is  very  large,  6F/F  -  -errcos  06(C  cos0/C-l). 

__  y 

With  our  assumption  that  Rq  is  centered  on  the  curve  C^cos0  =  C,  we 
have 


2 

6P  _  2e  cos  0 
? - 1 

Because  a  is  small,  this  ratio  can  easily  be  as  much  as  lOe. 
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IV. 


STATISTICAL  CONSIDERATIONS 


Thus  far ,  we  have  shown  how  to  calculate  the  backscattered  power 
averaged  over  a  statistical  ensemble  of  ocean  surfaces.  The  question 
of  just  what  "averaged  over  an  ensemble"  means  in  terms  of  physical 
measurements  remains.  This  will  be  our  next  topic. 


One  way  of  performing  an  average  is  to  look  repeatedly  at  dif¬ 
ferent  pieces  of  the  ocean  surface  which  are  far  enough  apart  to  be 
statistically  independent.  The  minimum  distance  between  two  such 
pieces  is  determined,  of  course,  by  the  correlation  length  of  those 
properties  of  the  ocean  surface  surface  which  are  important  in  the 
scattering  process.  Since  the  specular  part  of  the  scattering  cross- 
section  depends  only  on  the  mean  square  slope,  the  correlation  length 
relevant  for  specular  reflections  is  clearly  that  for  slopes,  which 
turns  out  to  be  some  tens  of  centimeters.  For  Bragg  scattering  (scat¬ 
tering  from  waves  of  a  definite  wavelength,  the  dominant  process  in 
backscattering  at  moderate  elevation  angles)  the  high-frequency  part, 
h^,  of  the  elevation  determines  the  scattering.  To  find  the  correla¬ 
tion  length  relevant  for  Bragg  scattering  we  must,  therefore,  study 
the  correlation  function 


P-^x-y) 


(h^xjh^y)) 


k-4ei£.(x-y) 


-4 

Because  it  contains  a  factor  k  ,  the  integral  on  the  right  is  rapidly 
convergent  and  receives  most  of  its  contribution  from  the  region 
kc  s  k£  2kc.  When  |  (x-y)  |  is  small  compared  to  Xc  =  2nk^,  the  ex¬ 
ponential  is  essentially  constant  over  this  region  and  p(x-y)  is  of 
order  2x10” However,  if  jx-y|  is  large  compared  to  Xc,  the  ex¬ 
ponential  factor  oscillates  rapidly  and  the  integral  is  very  small. 
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Thus,  we  may  take  the  correlation  length  for  Bragg  scattering  to  be 
something  like  which,  as  pointed  out  above,  will  be  a  few  radar 
wavelengths . 

Actually,  the  above  correlation  lengths  are  so  small  that  they 
are  of  little  interest  except  in  very  special  cases.  To  see  this  we 
have  to  understand  what  happens  when  a  radar  looks  at  the  ocean  surface. 

Suppose  a  radar  illuminates  an  area  A  of  the  ocean  surface.  Re¬ 
ferring  to  the  above  numbers,  we  see  that  in  general  A  will  be  very 
large  compared  to  the  relevant  coherence  length  squared.  Imagine  now, 
dividing  A  into  patches  whose  linear  dimension  is  of  the  order  of  a 
coherence  length.  We  can  write  the  backscattered  wave  as  where 
iL  is  that  part  of  the  backscattered  wave  which  comes  from  the  i-th 
patch.  Then  setting  ^  =  ai  ex^i ,  we  have 

p  _  a  a  ei(0£-0j) 

•  •  X  *1 

lO  J 

for  the  returned  power.  If  we  now  average  P  over  an  ensemble  of 
statisically  independent  areas  A,  the  averaged  power  is 

<r>  <a  a  eiCW>  <af> 

i,j  J  1 

which  follows  from  the  fact  that  the  phases  ex%  are  random. 

There  are  now  two  questions:  (i)  What  do  we  mean  by  independent 
areas?  and  (ii)  How  many  areas  are  needed  to  determine  (P>  to  a  given 
accuracy?  The  answer  to  the  first  question  is  almost  trivial.  In 
order  that  the  phases  e^i  be  uncorrelated,  the  two  areas  must  be  non¬ 
overlapping.  We  are  assuming  that  the  time  difference  between  mea¬ 
surements  is  less  than  the  decorrelation  time  of  the  phases.  To 
anser  the  second  question,  we  need  the  variance  of  P.  Here  we  appeal 
to  the  well-known  fact  that  for  a  sum  of  terms  wxth  random  phases, 
such  as  in  the  last  equation,  the  variance  is  always  of  the  same  o"der 
as  the  square  of  the  average,  i.e., 
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y<(p-<p»2>  =  <p> . 


This  means,  of  course,  that  co  measure  <P>  to,  say,  one  part  in  ten, 
we  need  one  hundred  separate  areas. 

Notice  that  the  above  conclusions  are  independent  of  the  size  of 
the  area  illuminated  by  the  radar  (so  long  as  its  linear  dimensions 
are  large  compared  to  the  correlation  length).  Thus,  contrary  to  one's 
first  impression,  the  accuracy  of  a  measurement  of  P  does  not  improve 
as  the  size  of  A  increases.  Also,  it  is  clear  that  the  magnitude  of 
the  small  correlation  lengths  does  not  enter  in  a  critical  ry. 

This  perhaps  surprising  situation  arises  because  a  radar  is  a 
coherent  source  of  radiation.  Suppose,  on  the  contrary,  that  the  source 
of  radiation  were  incoherent.  If  this  were  the  case,  the  equation 
given  on  the  previous  page  for  the  backscattered  power  should  be  re¬ 
placed  by 


p  =E 

i>j 


J 


where  0/  is  the  phase  of  the  incident  radiation,  assumed  to  vary  rap¬ 
idly  with  time  and  index  i  (this  is  to  incorporate  the  assumption  that 
the  incident  radiation  is  incoherent).  Averaging  over  a  time  long 
compared  to  the  coherence  time  of  the  0^  gives 


/ 

\ 


P) 


time 


average 


2 

The  point  is  now,  thav  the  average  of  E  a.  over  an  ensemble  of  areas 

*>  1  2 

A  gives  (P>  as  before,  but  the  vari  nee  in  E  a-  is  not  of  order  (P> 

2  i 

but  rather  of  order  (P)(L  /A)  ,  whtre  L  is  the  larger  of  the  coher¬ 

ence  length  of  the  radiation  and  the  relevant  coherence  length  of  the 
ocean.  Thus,  foi1  an  incoherent  source  the  accuracy  of  a  measurement 
does  increase  with  A. 
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The  reader  may  wonder  why,  in  the  case  of  a  coherent  source,  we 
did  not  average  the  power  equation  over  time  in  order  to  obtain  a 
result  similar  to  that  just  described.  The  reason  is  that  the  time 
scale  involved  is  vastly  different.  If  an  incoherent  source  has  a 
coherence  length  L,  its  coherence  time  is  very  small,  being  on  the 
order  of  L/c,  where  c  is  the  velocity  of  light.  With  a  coherent  source 
the  corresponding  time  is  the  coherence  time  of  the  ocean  surface. 

This  is  on  the  order  of  L'/v,  where  L7  is  a  coherence  length  for  the 
ocean  and  v  is  a  typical  wave  velocity.  It  is  the  large  ratio  c/v  — 

g 

10  that  makes  coherent  and  incoherent  sources  so  different. 

finally,  it  should  be  pointed  out  that  measurement  of  Bragg  scat¬ 
tering  has  some  statistical  properties  which  are  different  from  those 
of  specular  reflections.  Suppose,  for  example,  that  we  make  many 
measurements  of  specular  reflection  from  a  single  patch  of  ocean  sur¬ 
face  using  various  wavelengths  of  incident  radiation  and  angles  of 
incidence,  but  completing  all  the  measurements  within  one  coherence 
time.  The  statistics  of  the  measurement  have  not  been  improved  in 
this  case.  No  matter  what  angle  or  wavelength  we  use  to  measure 
specular  reflection  we  are  always  measuring  the  same  quantity,  namely 
the  mean  square  slope.  Thus  we  might  as  well  have  carried  out  all  the 
measurements  at  the  same  wavelength  and  angle,  gaining  no  improvement 
in  statistics.  Bragg  scattering  is  different,  however.  By  carrying 
out  the  measurement  at  different  angles  and  wavelengths  we  are  measur¬ 
ing  different  Fourier  coefficients  of  the  correlation  function  p^. 

Since  these  Fourier  components  are  statistically  independent  each 
measurement  gives  new  information  and  the  statistics  can  be  improved. 

As  an  example,  suppose  we  wish  to  measure  p^(0)  which  is  the  Fourier 
component  integrated  over  -C-space.  According  to  our  formulas  in  Sec¬ 
tion  III,  scattering  at  moderate  elevation  angles  samples  the  Fourier 

2  2  2 

components  of  over  an  area  of  order  ^racjar  m  in  t- space,  where  m 

is  the  mean-square  slope  of  the  ocean  surface.  Since  there  are  m-^ 

-2 

such  areas  available,  we  can  make  m  independent  measurements  whose 

o 

sum  (which  gives  p,  (0))  will  have  a  variance  of  m  times  the  variance 

A  A 

of  a  single  measurement.  Since  m^  is  of  order  10” ^  this  is  a  non¬ 
trivial  increase  in  statistical  accuracy. 
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FIGURE  2,  Partitioning  of  Scattering  Configuration 
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To  detect  ripples  on  the  surface  of  the  sea  by  passive  electro¬ 
magnetic  means  several  different  wavelengths  may  be  employed.  The 
question  naturally  arises  as  to  which  wavelength  range  is  most  suit¬ 
able.  The  answer  depends  not  only  on  the  physics  of  the  interaction 
of  electromagnetic  waves  with  the  surface  of  the  ocean  but  also  on 
the  technology  of  detectors.  In  this  section  a  comparison  of  the 
infrared  and  centimeter  wavelength  ranges  will  be  made. 

Several  physical  effects  contribute  to  the  detection  of  surface 
ripples  by  electromagnetic  waves.  Only  one  of  these  will  be  considered 
here.  A  detector  pointed  at  the  sea  surface  receives  reflected  radi¬ 
ation  from  different  portions  of  the  sky  due  to  the  presence  of  the 
ripple.  Since  the  radiance  of  the  sky  varies  with  elevation,  the 
presence  of  the  ripples  will  lead  to  an  average  variation  in  the  re¬ 
ceived  radiance  and  the  detection  of  the  ripples. 

Quantitatively,  the  spectral  radiance  (power  per  unit  area  per 
unit  solid  angle  per  unit  wavelength)  Wdgt  (\,  N^)  received  from  a 
direction  given  by  a  unit  vector  at  wavelength  \  consists  of  two 
parts:  (.1)  the  reflected  radiance  of  the  sky  at  a  direction  Ng  re¬ 
lated  to  by  the  law  of  reflection  and  (2)  the  emitted  radiance  of 
the  sea  itself.  In  terms  of  the  spectral  radiance  of  the  sea  W  , 
that  of  the  sky  W  ,  ,  the  reflectivity  p  of  sea  water,  the  zenith  di- 

SKy 

rection  z,  and  the  normal  to  the  sea  surface  N  we  write 


W 


=  P^-S-J!)Wsky  0„z-ns)  +  [l-pU,5.N>]WseaCX>  (1) 


The  angle  of  incidence  is  related  to  the  angle  of  reflection  by 


N 

~s 


(-1  +  N  N)N, 


(2) 
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where  I  is  the  unit  dyadic.  For  small  sea  slopes,  N  may  be  written 
as  z  +  s ,  Eq.  (1)  expanded  in  powers  of  £  and  averaged  over  the  dis¬ 
tribution  of  sea  slopes.  The  averages  of  e  and  may  be  expressed 
in  terms  of  the  r.m.s.  sea  slope,  m,  by 

<e>  =  -  m  z 


(3) 


Expressed  in  terms  of  the  zenith  angle  9  we  then  have  (see  Appendix) 
for  the  average  variation  in  received  radiance 

aw 


6Wdet(\,0)  =  m 


[ft 


) 


sky  dp 
+  "S0  W 


+  (W  .  -  W  ) 

*  sky  sea7 


( 


$  *  ”•  s 


)] 


(4) 


The  sky  radiance  contrast  arises  from  two  sources — absorption  and 
elastic  scattering.  Elastic  scattering  is  important  only  for  wave¬ 
lengths  £2u.  because  of  the  wavelength  dependence  of  the  elastic  cross 
section  and  the  size  distribution  of  the  scattering  particles.  In  that 
region,  the  sky  contrast  arises  because  on  the  average  the  reflected 
radiation  originates  one  mean  free  path  length  away  and  there  are  more 
scatterers  at  low  elevations  near  the  horizon  than  at  high  elevations 
near  the  zenith.  This  effect  can  lead  to  strong  radiance  contrasts  at 
low  elevation  angles  (see  Fig.  1)  of  the  order 


(scale  height  \  j 
for  scattererl  /  (mean  free  path) 
density  // 


^  1  km/10  km 


(Ip) 


Observations  at  such  small  angles  from  airplanes  are  different  because 
at  typical  airplane  heights  the  observation  distance  is  comparable 
with  the  mean  free  path.  We  will  not  consider  the  sky  constant  from 
elastic  scattering  further. 


no 


The  important  factor  contributing  to  the  sky  radiance  contrast 
is  absorption.  On  the  average  the  reflected  radiance  originates  one 
attenuation  distance  L(X)  away  at  a  height  L(X)  sin  (elevation  angle). 
Since  the  temperature  of  the  atmosphere  and  hence  the  radiance  varies 
with  height  this  will  lead  to  a  greater  radiance  at  the  horizon  than  at 
zenith.  The  resulting  contrast  will  be  small  for  those  wavelengths 
where  the  absorption  is  large  and  large  when  the  absorption  is  small. 

Figure  2  shows  the  experimental  sky  radiance  for  various  angles 
in  the  infrared  range.  No  contrast  is  observed  in  regions  of  strong 
absorption  (e.g.,  5(i  and  15(a)  while  the  maximum  contrast  is  obtained 
in  the  region  about  10(a.  In  figure  3  the  sky  radiance  (and  its  equiv¬ 
alent:  temperature)  from  this  data  at  10n  is  plotted  as  a  function  of 
angle.  Also  plotted  on  the  same  graph  is  the  spectral  radiance  for 
1.54  cm  normalized  to  the  same  height  at  6  =  90°.  Several  features  are 
clear.  Because  atmospheric  absorption  is  stronger  at  lOp  than  1.5  cm 
the  temperature  contrast  is  smaller  in  the  infrared  region  than  in  the 
microwave.  However,  because  the  dependence  of  radiance  on  temperature 
in  the  infrared  is  exponential  (hc/XkT  ~  5)  while  in  the  microwave  it 
is  linear  (hc/XkT  ~  1/50),  there  is  not  a  great  difference  in  the  ra¬ 
diance  constrasts. 

These  curves  are  the  first  elements  which  enter  into  a  calcula¬ 
tion  of  SW^  /Wd  The  second  element  is  the  reflectivity.  This  is 
estimated  from  the  standard  Fresnel  formulae.  For  the  X  =  1.5  cm  the 
curves  of  these  quantities  are  already  in  hand  in  Fig.  4 

The  largest  value  of  6Wdet  is  obtained  at  high  angles.  At  a 
typical  large  angle  of  8  =  75°  we  find  by  crudely  estimating  the  de¬ 
rivatives  of  these  curves 


det  _  .  2  .  ,.  .  .  _ 

— rj - w  0.6  m  ,  m  m  radians,  X  =  1.5  cm. 

oet 


Taking  a  temperature  resolution  of  0.2  K  we  have  approximately  for 
X  =  1,54  cm. 


Ill 


det  inn  2  .  .. 

-gTpy- —  «  100  m  ,  m  m  radians 

inst 

where  6W.  .  is  the  instrumental  radiance  resolution, 

inst 

For  the  10n  wavelength  the  index  of  refraction  was  computed  and 
plotted  in  Fig.  4  using  an  index  of  refraction  of  1.3.  The  curves  of 
W^et  for  horizontal  and  vertical  polarizations  are  given  in  Fig.  5. 
Crudely  estimating  F.q.  (4)  for  9  =  75°  one  finds 

6W 

-ii -e •  sa  0.6  m  ,  m  in  radians,  X  =  10p, 
aet 

The  similarity  of  this  number  with  that  obtained  for  X  =  1.54  cm  re¬ 
flects  the  similarity  of  the  sky  contrasts  at  the  two  wavelengths. 

If  we  take  6T  =  0.01°K  for  the  temperature  resolution  in  the 
infrared  we  find 


«  3000  rn 


2 


m  in  radians. 


The  conclusion  is  then  that  the  infrared  is  favored  over  micro- 
wave  radiometer  by  roughly  a  factor  of  30.  The  basic  reason  for  this 
is  that  the  window  at  10u  is  sufficiently  transparent  that  the  sky 
radiance  contrasts  are  almost  the  same  at  the  two  wavelengths,  while 
the  resolution  of  the  infrared  is  better  by  roughly  a  factor  of  30. 


112 


AFPENDIX 


DERIVATION  OF  EQUATION  (4) 

Let  n  ~  normal  to  sea  surface 

n  =  direction  to  observer 

n  =  direction  seen  on  sky 

~s 

Snell’s  law  is  expressed  by 

nxn  =  -  n  x  n 
~  ~o  ~  ~s 

Forming  the  vector  product  of  both  sides  with  n, 
n  x  (n  x  n  )  =  -  n  x  (n  x  n  ) 

r+f  /v  ~Q  ~  ^ 

whence 

n(n  •  n  )-n  =  -  n(n  -  nj  +  r. 

Since 


n  •  n  =  n  •  n 
~  “~o  ~  ~s 


we  have 


n  =  -  n  +  2n(n  •  n  )  =  (-  I  +  2nn)  •  n. 


rC*3  r^f  '  rsf  '■'-Q  '  '**'  '  ''"“O 

The  detected  radiance  is  equal  to  the  incident  radiance 

det  sky  ~  ~s' 

where  z  is  the  zenith  vector 

fV 

=  W  ,  (z  • (-1  +  2nn)  •  n) 
sky  ~  ~  o 

Now  for  small  sea  slopes  n  =  z  +  e, 

e  =  (sin0  coscp,  sin9  sincp,  cosQ-l) 
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or 


<e  =-  (Qcoscp,  Qsin^j  -  92/2)  +  0 CG 3 > 

Therefore,  accurate  to  second  order  in  9,  we  have 

Wdet  =  Wsky  *  C"~  +  22Z+  2(^-  +  +  2.S&]  * 

=  Wsky  *  So  +  2~  +  +  ^  *  So] 

=  W  ,  (z  •  n  )  +  W',  [2z  •  (ze  +  ez)  .  n  +  2z 

sky  ~  ~o  sky  ~  ~~  o  ~ 


+  2  W"  .  [z  .  (ze  +  ez)  •  n  3^  +  ... 
SKY  *** 


The  following  averages  will  be  needed: 

<£>  -  -  \  <e2>  z 


(ee)  •  n  1 

r'^r+J*  (vQ 


< ££>  =  al  +  b  zz 
a  =  ^  <92) 

a  +  b  =  0  b  =  -  k  <92> 


Using  these  results,  the  first  bracket  in  the  expansion  of  wdet  becomes 
on  averaging  -2  <02>  (z  •  n  ).  The  second  bracket  is  equivalent  to 


2[e(I  +  zz)  .  n^]2  =  2^  .  (I  +  zz)  •  ee  .  (I  +  zz)  . 

=  <e2>  n,  *  (I  +  z&)  *  (I  ~  zz.)  *  <3  +  J5S)  *  2, 

=  <&2>  So  ‘  CJ  -  z&)  •  no 

*  <e2>  [1  -  (nQ  •  z)2]  =  <02>  sin29o 
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where  9  is  the  angle  of  observation  measured  from  the  vertical, 
o 

^det  =  ^(cos®0)  +  (Q2)  C-2W'  cos0o  +  W"  sin2eQ] 


Now 


W'  = 


dW  1 
d0  sinff 


W"  =  + 


_ 1_  d_  /dW _ 1 

sin9  d9  IdF  sxn 


sm$J 


1  cose  dW 

^  de2  “  sin3e  m 


Thus, 


w  =  w(9  )  +  <e2)  +  cote  3W 

det  o  gg z 


This  computation  was  made  leaving  out  the  reflectivity.  Actually: 


Wdet  =  ^ 


iV  sky 


~s 


Expanding  the  reflectivity  in  the  same  manner  as  before, 


p(n  •  n  )  =  p(z  •  n  +  e  •  n  ) 
~o  ~  o  ^  ~o 


=  p(z  .  nQ)  +  p'(e  •  nQ)  +  \  p"(e  •  nQ)2  +  ... 

Including  the  extra  terns  in  the  expression  for  t 

W  =  [p  +  p '  ( e  •  n  )  +  (2g)  pw(e  •  n  )2]  |  W  , 

det  y  ~  o  2  v  ~  ~o  I  sky 

+  W',  [2z  .  (ze  +  ez)  •  n  +  2z  (ee)  n  1 

+  2Wsky  *  nJ2! 


Thus, 
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The  net  result  is  that  the  variation  in  radiance  is  given  by 


This  expression  accounts  for  the  contribution  to  the  variation  in 
received  radiance  due  to  the  first  term  in  Eq.  (1).  The  contribution 
of  the  second  term  can  be  obtained  by  inspection  of  the  above  expres¬ 
sion  (replacing  WsRy  by  Wgea  and  noting  that  Wsea  is  independent  of  0) 
Thus,  the  contribution  of  the  second  term  is 


6W 


det 


-(m2/4)  W 


sea 


.2 

— £  +  cot9 

ae 


dp 

TO 


Combining  these  results  yields  Eq.  (4)  of  the  test. 
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FIGURE  1 .  Spectral  radiance  of  a  clear  sky,  showing  the  dependence  of  the 
scattered  radiation  on  elevation  angle.  Measurements^made  from  Colorado  Springs, 
Colorado,  near  noon  in  September;  elevation  angles  0  (top  curve),  7.2  and  30 
(lowest  curve).  Note  that  the  ordinate  scale  for  the  short-wavelength  set  of  curves 
is  lOtimes  larger  than  the  scale  of  the  longer-wavelength  curves.  (From  Bell  et  al 
0960] 
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FIGURE  2.  Spectral  radiance  of  a  clear  sky  (for  angles  of  elevation  0,1.8, 
3.6°,  7.2°,  14.5°,  30°,  and  90°  above  the  horizon.)  Ambient  temperature  8°C; 
measured  at  night  in  September  from  Elk  Park  Station,  Colorado,  11,000  ft  above 
sea  level .  Ambient  temperature  about  2/  C;  measured  in  June  from  Cocoa  Beach, 
Florida.  (From  Beil  et  al  0960] 
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REFLECTIVITY 


FIGURE  4.  Reflectivity 
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Sophisticated  signal  processing  is  claimed  by  some  of  its  advo¬ 
cates  to  make  possible  remarkable  improvements  in  detection  systems. 
While  this  is  true  in  some  cases,  in  others  the  claims  are  extravagant. 
Tlris  note  is  written  to  provide  a  basic  criterion  for  detection  that 
depends  only  on  hardware  capability  and  by  which  ultimate  performance 
limits  can  be  set:  limits  that  can  be  approached  but  not  surpassed 
by  astute  processing  and  presentation.  The  basis  for  this  is  the  work 
of  Harris*  on  decision  theory.  The  derivation  of  one  of  his  main  re¬ 
sults  is  abstracted  here  for  completeness. 

We  will  look  at  the  detection  problem  as  a  binary  decision  between 
two  signal  sources  denoted  I  and  II,  with  the  signals  accompanied  by 
additive  Gaussian  noise.  For  definiteness,  take  a  two-dimensional  data 
presentation  with  mean  flux  densities  of  H^^y)  and  Hjj(x,y)  for  the 
sources  I  and  II.  The  likelihood  that  a  set  of  flux  readings  R, ,  R,. 
...R  with  dispersions  cr^,  ...  cT  for  patches  of  area  AxAi?  is  ob¬ 
served  in  response  to  source  I  is 


n 


L(I)  =  7 r 


expt-CRj  -  Hn  Ax4y)2/2oJ] 


where  R.^  is  the  flux  measured  at  the  display  point  (xi,yi),  and  HIi  = 
Hi(xi,7i).  Similarly,  the  likelihood  of  the  same  readings  in  response 
to  source  II  is 


L(II)  =  7T  - 

i=i  y[: 


exp[- 


2 

2TT  a 7 
l 


(Ri  - 


H. 


Hi 


AxAy)2/2c?] 


This  formulation  assumes  that  the  incremental  area  AxAy  is  sufficiently 
large  for  the  observed  readings  R^  to  be  regarded  as  statistically 
independent . 
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For  the  purpose  of  this  discussion,  the  relative  risks  of  error 
in  deciding  that  soui'ce  I  is  present  when  source  II  is  really  present, 
and  vice  versa,  can  be  ignored.  Under  these  conditions,  the  decision 
rule  is  to  select  the  alternative  with  the  larger  likelihood.  In  terms 
of 


t  =  21og[L(I)/L(II)] 


the  procedure  is  to  decide  that  I  is  present  if  i|r  > 
that  II  is  present  if  •)  <0.  Specifically, 


=  £  u/opt-a^Hjn 

l=x 


HIjL)  AxAy  +  (Hj. 


Now  suppose  that  source  I  is  actually  present; 
R.  =  HIi  Ax  Ay  +  ni 


0,  and  to  decide 


HIJ)(Ax)2(Ay)2] 

then 


with 


c2 

1 


v^AxAy 


Here,  m  represents  the  additive  Gaussian  noise,  and  v^  is  the 
noise  variance  per  unit  area.  Substituting  this  into  the  expression 
for  *  yields 


<Hii  - 


v. 

i 


f  2ni  (Hiii  -  hi? 
i=l  vi 


The  mean  of  Vj  can  be  expressed  as 


=// 


[HTT(x,y)  -  HT(x,y) 1' 


1 Illlil  "  I 

v(x,y) 


dxdy 
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The  integral  is  taken  over  the  area  of  presentation  (e.g.,  the  field 
of  view  of  display).  The  variance  of  4'j  is  calculated  to  be 


The  probability  of  correct  decision  (i.e.,  that  >0)  is 


w 

/2 

e-2  ^  dz 


(1) 


-U./E 


The  probability  of  a  correct  decision  depends  on  the  single  parameter 


H/2  = 


(%) 


(2) 


when  N  denotes  the  background  noise  per  unit  area. 

Well-matched  processing  and  data  presentation  can  take  full  ad¬ 
vantage  of  the  signals  provided  by  detection  equipment,  but  cannot 
improve  the  probability  of  a  correct  decision  over  that  implied  by 
Eq.  (1)  and  (2). 
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ii  abstract 


(S^This  study  deals  with  fundamentals  in  the  performance  of  airborne  sensors  for  detecting  the  wake  frcra 
the  oassagc  of  a  sufrarlne  through  stratified  water.  No  attempt  has  beer,  made  to  inbed  the  theoretical  results 

ir  .m  ope  rat  i  or.  a  ^  setti  g^  ge*^ated  by  fluid  displacement  idipole  source)  and  wake  collapse  (quadruple  source) 
The  wale  source  predominates  far  astern  frcn  shallow  moving  bodies,  an..  vice  versa.  The  associated  sur.ace 
straits  lead  to  a  characteristic  pattern  of  enhancement  of  ambient  wind  waves,  and  this  ir.  turn 


natural 

Kw-i  rfr<J  tnstrtsrent  noise  level,  the  detection  of  an  anomalous  pattern  covering  10m  x  1000m,  say,  re- 

gui  re '.strains  of  5(10-0  for  some  passive  (X?  and  micrcwave)  and  active  (radar)  devices  (but  the  limitations 
for  three  methods  differ  distinctly).  tv  , 

i)  There  is  soro  indication  tnat  the  detectable  strain  might  be  reduced  tc  0(10“5)  for  an  active  device 
0«  .  •),  a-ulai  resol utvor  ar.d  frequency  (Lasei)  viewing  the  sea  surface  2  ros  slopes  from  the  vertical.  Active 
dev:  as  appear  attractive  also  for  ere  reason  that  they  eliairate  the  last  step  in  an  already  most  involved 
seric-  of  c  >  .pied  phenomena:  internal  uave  generation  -  surface  straining  -  differential  roughe  -rg  (affec.irg 
scattering)  -  radiorc-nc  temperature  (IK). 
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